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INFELD FACTORIZATION AND ANGULAR MOMENTUM! 


By H. R. Cots 


ABSTRACT 


The connection between Infeld factorization operators and angular momentum 
operators, well known for spherical harmonics, is extended to other factorization 
problems by explicitly recognizing them as angular momentum problems. These 
other problems are: the symmetric top, electron-magnetic pole system, Wey]'s 
spherical harmonics with spin, free particle on a hypersphere. The Kepler 
problem is also included for it may be thrown into the form of a four-dimensional 
angular momentum problem. The transformation to momentum space for this 
problem is very much simplified by the connection between Infeld factorization 
and angular momentum. 


1. INTRODUCTION 


There are two main papers on the Infeld factorization method (Infeld 1941; 
Infeld and Hull 1951). The present paper partially answers the question 
raised by the suggestion made by Infeld (1941) that the factorization method 
is something more than merely a mathematical trick. 

In the problem of the spherical harmonics,* it is clear that the Infeld 
factorization operators are closely related to the angular momentum ladder 


operators, 
ka A = 
(1) L,+1L, =e Yi cot 9 as = 6, ‘ 
With P,"(cos #)e** written for the operand this ladder property yields 
J ‘ d \ m m+1 
(2) \™ cot 6 + aos P,"~P, 


in agreement with the recurrence relations obtained by the Infeld method. 
Similarly, by elimination of the co-ordinate ¢ from the angular momentum 
identities 
(L,+1L,)(L,—1L,) = L’-L,(L,+1), 


(3) AHA 0 ODD, 


the doubled factoring characteristic of the factorization method is obtained. 


1Manuscript received October 7, 1955. : ; mv d 
Contribution from the Department of Mathematical Physics, University of Manitoba, 


Winnipeg, Manitoba, and the Canadian Mathematics Congress Summer Research Institute, 
Queen's University, Kingston, Ontario. 

*Infeld and Hull (1951) Sec. 4.1. 

tInfeld and Hull (1951) page 30, footnote 6. 
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2. THE SYMMETRIC TOP,* THE ELECTRON-MAGNETIC POLE,f WEYL’S 
SPHERICAL HARMONICS WITH SPINt 

These are problems for which the components of the total angular momen- 
tum are constants of the motion. Although the connection has not before been 
pointed out explicitly, one might expect that again the angular momentum 
ladder operators would lead directly to Infeld factorization operators. 

In the symmetric top problem we use Euler angles such that 6, ¢ are the 
polar co-ordinates of the axis of symmetry about which the spin y takes place. 
We find for the total angular momentum 
(4) J =L+a 
where 

az = py csc A cos ¢, 


(5) a, = py csc Osin ¢, 
a, = 0, 


py being the momentum canonically conjugate to y, while the components of 

L have exactly the form of the components of orbital angular momentum 

of a particle and hence as operators will have the usual orbital angular momen- 
tum operator form. 

The ladder operators are 

: JrriJ, = L,t1Ly+a;+iay. 


Operating on an eigenfunction of L, = —1 0/0@ for the eigenvalue m these 
take the forms 

= sip Sn lu ae __ Pe = 

(6) J,+1J, e \" cot 6 a + 36 


which give the Infeld factorization operators for this problem. 
From the identities corresponding to (3) one could get the Infeld double 





factoring. 

The energy eigenvalues in terms of the moments of inertia A, C 
; oe FGAY = 
(7) ee ae 


follow immediately from the eigenvalues J(J+1) for total angular momentum 


and K for py. 

Similarly the Infeld factorization may be derived from angular momentum 
properties and the ladder operators related to angular momentum operators 
in the electron-magnetic pole problem and in the case of Weyl’s spherical 
harmonics with spin, if in the former the electromagnetic angular momentum 
and in the latter the Pauli spin are taken into account. 


3. FREE PARTICLE ON A HYPERSPHERE § 


The problem of the motion of a free particle on a hypersphere leads to 
Infeld’s first example of the factorization method. It is an angular momentum 


*Infeld and Hull (1951) Sec. 4.5. 
tinfeld and Hull (1951) Sec. 4.7 and Dirac (1931). 
tInfeld and Hull (1951) Sec. 4.6 and Weyl (1931). 
§Infeld (1941) Sec. 1 and Sec. 7. 
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problem but in four dimensions, and this occasions a change in the style of 
argument. There are now ‘wo angular momentum vector operators, namely 


(8) =-irxXV, 
J 0 
9 = - ae! 
(9) a) x4V —) 
where the vector notation denotes vectors under rotations of the x1, x2, x3 
space. 


L is just the usual three-dimensional angular momentum vector operator 
satisfying the relation 
(10) LXL=iL. 
Other angular momentum relations which may be verified by using (8) and 
(9) are 
(11) DxD=ilL, 
(12) DXL+LxD=2iD. 
These relations hold for a system of particles as well as for a single particle, 
but a fourth relation exists for a particle which does not hold generally for a 
system of particles 
(13) D-L=L-D=0. 
In terms of these operators the Schrédinger equation on the hypersphere of 
radius R 


(14) Ay+2mEy = 0 
becomes 

(15) (D?+L?) ¥ = Ay 
where 

(16) \ = 2mE R’. 


Since L? commutes with D? we may impose 
(17) L’y = l(l+1)y¥ 
and our problem becomes one of finding four-dimensiona] angular momentum 


ladder operators which raise and lower /. 
To find these and relate them to Infeld factorization operators we introduce 


the Pauli spin matrices o1, 2, ¢3, and the vector 


(18) 6 =iotjot+ko; 
with the properties 

(19) 6X 6 = 214, 

(20) (é-a)(é-b) = a-b+i6-a X b, 


where a, b are two vector operators which need not commute with each other 
but must commute with 4. 

The operator 
(21) K =6-L+1 
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has the property that 
(22) (K-1) K =L? 


which follows from (10) and (20), so that an eigenfunction of K for the eigen- 
value k is an eigenfunction of L? for the eigenvalue /(/1+1) with / = k—1 if 
k is positive and / = —k if k is negative. We take y to be a two-component 
column matrix which is an eigenfunction of K so that both its components are 
eigenfunctions of L? for the same value of /. 

Any operator which anticommutes with K changes y into another eigen- 
function of K for the eigenvalue K = —& and hence into another eigenfunction 
of L? for the eigenvalue (/+1)(/+2) if k is positive and for the eigenvalue 
(/—1)l if k is negative. That is, it raises or lowers / by one depending on 
whether & is positive or negative. It is a ladder operator. 

Just this property of anticommuting with K is held by the operator 
é-D as follows from (12), (13), (20). Furthermore as a consequence of (11), 
(20), (22) we have 
(23) D?+L? = (6-D)?+K?-1 


so that 6-D commutes with D?+L?. Hence it is exactly the ladder operator we 


seek. 
In terms of polar co-ordinates a, 6, ¢, given by 


x; = Rsina sin 6 cos ¢, 


x2 = Rsina sin 6 sin ¢, 
(24) ; 
x3 = Rsinacos 8, 
x, = Rcosa, 


the ladder operator é6-D becomes 


(25) é-D= =te4 2 (K—1)cot «| 
where 

(26) or = (6-r)/r. 
Applying this to a y taken to be of the form 

(27) vi = Yi, ) ur(a) 


where Y; is a (matrix) eigenfunction of K for k = /+1 whereas wu, is an ordinary 
function, we obtain 


(28) ~iey| £1 cot J] uj ~ Vig tus 


Cancelling off the factors depending on 6, ¢ only we have 
(29) = cot a Uy™ Uj41- 
Similarly by taking Y; to be an eigenfunction of K for k = —/ we get 


(30) [s+ (1+1) cot al U;,~ Uy}-1. 
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These are just the recursion relations which would be obtained by the Inield 
factorization method. 
We can derive the doubled factoring of the Infeld method by considering: 


(31) (6-D)? = D?+L?— K?2+1 


or 
(32) =a aw cota | o| 2—(K-1) “ni “ sh [axis |v. 
Oa 0a 


Now we permute the oa, in the middle through to the left to combine it with 
the o, already there. Making use of the facts that o, commutes with 0/da 
and anticommutes with K, and that its square is unity, we obtain: 


(33) | (x41) cot eet (K—1) ea y= [ a—xe-+1 |v. 


On substituting an eigenfunction of AK for k = /+1 and cancelling off the 
6, @ dependent factors the result is 


(34) | +2) cot a+ 4], cot rt i; = [ a—10+2) | Uy. 
da da 
On the other hand the use of an eigenfunction of K for k = —/ yields 
(35) eS pei [ /-+1) cot a+ : ] 4, = | a a-no+n | Uj. 
da da 


These last two equations form just the factorings of the factorization method 
and here have been obtained by angular momentum arguments. The factori- 
zation operators have been related to four-dimensional angular momentum 
ladder operators. 

Similarly the Infeld factorization for the problem of the particle in a space 
of constant negative curvature could be related to angular momentum 
properties. 

The factorization of the Bessel equation may be obtained from the above by 
taking the limit of zero curvature. 

It must be noted that the argument in this four-dimensional case depends 
essentially on a relation (13) which holds for a single particle but not in 
general for a system of particles. Hence it would not be possible to extend this 
treatment as can be done in the three-dimensional case for the symmetric 
top. 

4. THE KEPLER PROBLEM 


Jackson (1953) has pointed out that the Kepler problem can be thrown into 
the form of a four-dimensional angular momentum problem with 


(36) D = (—2mE)[}(p X L — L X p) —me’r/r), 
(37) 2E(D?+L?+1) = —me’, 


where D and L have the four-dimensional angular momentum properties 


(10), (11), (12), (13). 
Now the ladder operator from the angular momentum point of view is 
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(38) 6-D = —(—2mE)-[(6-V)K+me’s,], 

where we have made use of 

(39) (6-V)K—K(6-V) =ie-(V X L—L XV) 

from (20) and (21), and also of the fact that 6-V and K anticommute. Putting 
our ladder operator in terms of polar co-ordinates we finally obtain the result 


(40) ro (—2mBy 4%] (2-£=1)e+me'| 


We next apply this to an eigenfunction of A for the eigenvalue k = /+1, 
that is to an operand of the form (27) but with u;(a) replaced by a function of 
the radial co-ordinate, 
ee | 0 i 2 
(41) 6-D = —(—2mE) ‘oc, as (/-+1)+me* |. 
From this we can pick out the operator on the radial function which is clearly 
just of the Infeld type for raising / 
(42) 1 
dr r 1+1° 

Similarly we may pick out the Infeld ladder operator for lowering / by choosing 
k= -l 

d ,l+l1 me” 
(43) _ r b 
Thus the Infeld factorization of the Kepler problem is also a consequence of 
angular momentum properties. 


5. KEPLER PROBLEM IN MOMENTUM SPACE 


A problem to which considerable effort was devoted in the thirties is that of 
finding the momentum space wave functions for the Kepler problem (Podolsky 
and Pauling 1929; Hylleraas 1932; Elsasser 1933; Fock 1935). With the 
ladder operators as we now have them in terms of physical operators it is 
possible to obtain the result in a few lines with elementary operations. 

The hamiltonian for this problem is 


(44) (p?/2m)—e?/r = E 

so that in momentum space the operator e*/r may be represented by 
(45) e?/r = (p?/2m)—E. 

If we write 

(46) eo, = e'6-r/r = 1(6-V,)(p?—2mE) /2m, 


where V, is the gradient operator in momentum space, we may transform 
é-D (38) into a momentum space operator, 
(47) 6-D = —i(—2mE)-}[(6-p)K +(6-V,)(p?—2mE)/2]. 


Introducing polar co-ordinates p, 6), ¢p, in momentum space we may write 
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0 K-1 
i wna Lap “ae 
where 
(49) op = (6-p)/p. 


We now substitute this into (47) and pick out the recursion relations for the 
‘radial’ function P(p) in momentum space, 





P d apd+1) = 
(50a) ' pt p—2ml InmE (p° —2mE) Pi~ Put 
, I+1__ 2pl = _ 


The case / = n—1 is the ‘top of the ladder’ and for this case the operator of 
(50a) annihilates the corresponding P, 


s d_n-l, 2np_ | 
(51) E ? p—2mE (p” —2mE) P,-1 = 0. 


This is a first-order differential equation for P,-1 whose solution is 


n—1 
(52) Pra ~ (p> —2mE)"* 
The solutions for other values of / may be obtained simply by applying the 
operator of (50) for lowering /. 
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HIGH RESOLUTION RAMAN SPECTROSCOPY OF GASES 
VI. ROTATIONAL SPECTRUM OF SYMMETRIC BENZENE-d;' 


By A. LANGSETH? AND B. P. STOICHEFF 


ABSTRACT 


The pure rotational Raman spectrum of Ces6HsD3 vapor at a pressure of 
15 cm. Hg was photographed in the second order of a 21 ft. grating. The value 
of the rotational constant was found to be By = 0.1716;+0.0001 cm. This result 
confirms the earlier spectroscopic values of the internuclear distances in the 
benzene molecule. 


A. INTRODUCTION 

Recently the carbon-carbon distance in the benzene molecule was evalu- 
ated by two methods, each capable of high precision, but the resulting values 
were significantly different. From the rotational Raman spectra of gaseous 
CeHs and CsDe, Stoicheff (1954b) obtained the value of 1.3973;+0.001 A for 
the carbon-carbon distance in benzene while from the X-ray diffraction 
analysis of crystalline benzene at —3° C., Cox and Smith (1954) obtained the 
value of 1.378+0.0033 A for this distance. These values differ by 0.02 A, that 
is, by approximately ten times the quoted errors. More recently, Cox, Cruick- 
shank, and Smith (1955) have been able to ascribe most of this difference to a 
systematic error in the X-ray analysis arising from the anisotropy of the 
thermal oscillations of the benzene molecule as a whole. The corrected X-ray 
value of the carbon-carbon bond length is 1.392 A which is now in satisfactory 
agreement with the spectroscopic value. 

Even though there is agreement at present between the X-ray and spectro- 
scopic values, it seemed worth while to continue with an experiment conceived 
prior to the explanation of the original discrepancy, in order to check the 
spectroscopic value. The present investigation of the rotational Raman 
spectrum of symmetric benzene-d; confirms the earlier spectroscopic results 
obtained with CgHg and C,Ds, namely that the internuclear distances in the 
benzene molecule are 


ro(C-C) = 1.397; A and ro(C-H) = 1.084 A. 
B. OBSERVED SPECTRUM 


The apparatus has already been described (Stoicheff 1954a). A cylindrical 
lens was mounted in front of the photographic plate resulting in a 10-fold 
reduction in exposure time with no noticeable loss in resolving power. The 
sample of CsH;D; had a purity of about 90 mol. %, the remainder being 
mainly 1,3-CsH,D» with a trace of CsHsD. The pressure of benzene was 15 cm. 
Hg and the exposure time one to two hours. 

The rotational Raman spectrum of CsH;D; similar to that of CsH¢ consists 
of a series of sharp lines extending to about 50 cm on both sides of the 


1Manuscript received December 14, 1955. 

Contribution from the Division of Pure Physics, National Research Council, Ottawa, 
Canada. Issued as N.R.C. No. 3874. 

*Chemistry Department, University of Copenhagen; visiting Professor at N.R.C. in 1954. 
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exciting line. These are the lines of the S branches for which AJ = 2. The line 
spacing is therefore 4B, and is observed to be about 0.68 cm~'. The lines of 
the R branches (AJ = 1) have half the spacing of the main series of lines and 
are unresolved, thus producing a continuum which extends about 20 cm~! 
from the exciting line. 

Two plates of the CsH;D; spectrum were evaluated. About 120 Stokes and 
anti-Stokes lines were measured on each plate and the agreement between the 
line measurements on the two plates was within +0.02 cm—'. The average 
values of the wave number shifts are given in Table I. 


TABLE I 
OBSERVED WAVE NUMBER SHIFTS OF THE ROTATIONAL RAMAN LINES OF CeH;3D3 











Ff | Av |,2 c—o, FF | Av |,@ c—0, 

cm! em! | cm7 cm™ 
21 15.44,° +.003 51 36.035 ~.003 
22 16.114° +.019 52 36.72 —.007 
23 16.812 +.008 | 53 37.435" —.032 
24 17.506 +.003 | 54 38.07,* +.017 
25 18.20;* —.008 | 55 38.751 +.023 
26 18.87 6*° +.003 | 56 39.455 +.002 
27 19.54;*° +.019 57 40.132 +.013 
28 20.24, +.005 58 40.835 —.004 
29 20.93: +.001 | 59 41.51: +.004 
30 21.63; —.006 | 60 42.196 +.006 
31 22.305 +.006 61 42.875 +.009 
32 23.000 —.003 62 43.575 —.006 
33 23.685 —.005 | 63 44.26; —.010 
34 24.3% —.024 | 64 44.94, +.001 
35 25.062 —.006 | 65 45.636 —.007 
36 25.747 —.005 66 46.335 —.021 
37 26.445 —.012 67 47.02; —.023 
38 27.125 —.011 | 68 47.679* +.004 
39 27.805 —.004 69 48.357" +.012 
40 28.47;* +.014 70 49.065 —.006 
41 29.122* +.050 | 71 49.749 —.001 
42 29.849 +.018 I 72 50.446 —.016 
43 30.54 +.002 | 73 51.126 — 017 
44 31.215 +.016 74 51.79; —.003 
45 31.905 +.011 i 75 52.48> —.004 
46 32.59. +.010 | 76 53.138 +.025 
47 33.27; +.012 | 77 53.833 +.015 
48 33.972 +.002 i 78 54.516° +.016 
49 34.66, —.005 | 79 55.227° —.010 
50 35.35: —.006 | 80 55.903° —.001 





2A]] values are averages of Stokes and anti-Stokes lines measured on two plates, 
except the values marked °. 

Measurement on one plate only. 

*Lines blended with grating ghosts. 


C. ROTATIONAL ANALYSIS 


The method of analysis used to determine the rotational constants is de- 
scribed in section D of paper II (Stoicheff 19546). Once the rotational number- 
ing was established, a graph of the values | Av |/(J+3/2) was plotted against 
(J+3/2)?. The intercept on the ordinate axis yielded a value of 4By and the 
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slope a value of 8D,. From these values, the rotational constants 
By = 0.1716;+0.0001 cm-! and Dy, = 1.3X10 cm“! 


were obtained for the CsH;D; molecule. As a check, these values were used 
to calculate the wave number shifts of the Raman lines which were then 
compared with the observed shifts. The agreement between the calculated 
and observed values is within +0.03 cm! for all the lines, as shown in Table I. 
The effective moment of inertia about any axis in the plane of the molecule 
passing through its center was calculated from the above value of By and the 
atomic constants of DuMond and Cohen (1953), giving 
T;° = (163.06+0.10) K10-*° gm. cm? 


Within the experimental error, this value is equal to the mean of the moments 
of inertia of CsHe and CsDs, 147.59X10-*° and 178.45X10-* gm. cm.?, re- 
spectively, thus confirming the earlier results obtained from the Raman 
spectra of CsHe and CsDsz. 

The agreement between the results of the present investigation and the 
earlier one can also be seen from a graphical determination of the internuclear 
distances. In Fig. 1 are shown the graphs of the C-H distances and C-C dis- 


1.07 





Fic. 1. Graphical determination of the internuclear distances in the benzene molecule. 
The graphs represent the values of the C-H and C-C distances that are consistent with the 
measured values of Bo for the molecules CeH., CoD, and CsH3D3. 


tances consistent with the By values of CsHe, CsDs, and CsH3;D3, assuming a 
planar hexagonal structure. The values of the internuclear distances are given 
by the intersections of any two of the curves. All three sets of values are in 
close agreement. 

A least squares calculation using the three moments of inertia to determine 
the two internuclear distances of the benzene molecule gives the values 


ro(C-C) = 1.397,40.001 A and 79(C-H) = 1.084+0.005 A. 
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Although it is fortuitous that these values agree exactly with the earlier values, 
it seems safe to conclude that the spectroscopic values are correct to within 
the quoted errors. 
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ON THE THEORY OF A COAXIAL TRANSMISSION LINE 
CONSISTING OF ELLIPTIC CONDUCTORS! 


By J. Y. WonG 
ABSTRACT 


A coaxial transmission line consisting of conductors of elliptic cross-section 
is treated as a boundary-value problem. Using elliptic-cylinder wave functions, 
a formula is derived for the propagation constant of the allowed symmetric modes 
of propagation. A condition for the generation of the higher-order modes is 
established in terms of the physical parameters. The analysis is restricted to the 
case where both the inner and outer conductors are confocal. A configuration of 
practical interest follows for the limiting case of the inner conductor consisting 
of a flat strip. The resulting structure can be regarded as constituting a form of 
shielded “‘strip’’ transmission line. The analysis may be used to provide an 
approximate theory for the ‘rectangular’ coaxial line. 


INTRODUCTION 


In this paper, a coaxial transmission line consisting of conductors of elliptic 
cross-section is treated as a boundary-value problem. Using elliptic-cylinder 
wave functions, a formula is derived for the propagation constant of the allowed 
symmetric modes in the direction of propagation. The analysis is restricted 
to the case where the inner conductor and the outer conductor are confocal. 
An expression for the critical frequencies of the higher-order modes is found 
to be expressed in terms of a simple function of the physical parameters of 
the line. An obvious extension of the present analysis is the determination of 
the characteristic impedance. 

A configuration which may be of practical interest results for the limiting 
case of the inner conductor consisting of a flat strip. For large ellipticities of 
the outer conductor, the structure may be regarded as constituting a form of 
shielded ‘‘strip’’ transmission line. A common type of strip line encountered 
in practice consists of a flat conductor placed between two parallel conducting 
ground planes in the form of a sandwich. Papers describing this type of line 
have appeared extensively in the literature (Begovich and Margolin 1950; 
Cohn 1954; Trans. I.R.E. 1955). It is apparent that such a configuration 
cannot be analyzed rigorously from the electromagnetic point of view and the 
determination of some of the electrical parameters must be restricted necessarily 
to a study of the electrostatic model. In the present case the electric field is 
confined entirely within the region between the conductors; consequently 
loss due to radiation and a reduction in the characteristic impedance due to a 
fringing field do not exist, as in the case of the conventional shielded strip 
line. Finally, the present analysis may be used to provide an approximate 
theory for the rectangular coaxial line. 


WAVE EQUATION IN ELLIPTIC CYLINDER COORDINATES 


In problems dealing with cylinders of elliptic cross-section, it is necessary 
to introduce an elliptic cylinder coordinate system. Elliptic cylinder coor- 
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Contribution from the Radio and Electrical Engineering Division, National Research 
Council, Ottawa, Canada. Issued as N.R.C. No. 3878. 


354 








WONG: COAXIAL TRANSMISSION LINE 355 


dinates (u,v, z) are related to cartesian coordinates (x, y, z) through the follow- 
ing transformation equations: 


x = Cy cosh u cos 2, 


(1) y = cosinh usin v. 


The coordinate surfaces are confocal elliptic and hyperbolic cylinders, 
2 2 


x y 

cotcosh's * cosinh's wn 

(2) we 

cocosv cosin'y 

and planes z = constant. From eq. (2) it is apparent that co is the semifocal 
distance of the ellipse. 

Fig. 1 illustrates a family of confocal elliptic and hyperbolic cylinders. 








eT 
COTTE S 






Fic. 1. Elliptic cylinder coordinate system. 


In a cylindrical coordinate system, the electromagnetic field can be resolved 
into two partial fields, each derivable from a purely scalar function II, and 
II,* respectively. Let y represent either function, then y must satisfy the wave 


equation 

(3) (V?+k?)y = 0, 
where V? is the Laplacian operator, and 

(4) k? = —iwpu(otiwe). 


Here the rationalized m.k.s. system of units is employed and the time function 
e“ is understood. In eq. (4), u is the permeability, ¢ is the conductivity, and 


e is the dielectric constant of the medium. 
In elliptic cylinder coordinates, elementary solutions of eq. (3) can be 


represented by the form 


v = f(u,v) et”, 
where f(u, v) is a solution of 


(5) a5 42 OS ea? I’ (cosh’u—cos yf = 0; 


in eq. (5) 
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(6) r? = k?—h?. 


Using the method of separation of variables, eq. (5) yields the two well-known 
radial and angular Mathieu equations. Thus, within a homogeneous isotropic 
domain, every electromagnetic field can be represented by a linear combination 
of the elementary wave functions 


~ ih 
) We, Om = Je,0m(s, U)S€,0m(s, v) eX, 


~ 
‘ 


( 
(8) We, On = 118.’om(s, U)S,0m(S, v) et. 


In eqs. (7) and (8), s = co?T*. The wave functions consist of both even and 
odd cylinder functions and, of these two equations, (7) applies to the finite 
domain. At great distances from the source eq. (8) must be employed, since 
it reduces asymptotically to a wave travelling radially outward. 

Definitions of the various radial and angular wave functions can be found 
in numerous sources in the literature (Morse and Feshbach 1953; Sinclair 
1951; Stratton, Morse, Chu, and Hutner 1941), and will not be repeated 
here. 

Consider an elliptic cylinder u = uo of infinite length, immersed in an 
infinite homogeneous medium. Utilizing the results of eq. (7), the field com- 
ponents at all interior points, « < uo, can be written down as follows: 


co 


(9) Ey, = = ~ (ih/T17hy) Je,0m(S1, U) S€,0m(S1, ¥)Ae,Om 
m=0 
— (iw /T17h,) Je,0m(S1, U) S€,0m(S1, 0)b€,0m} e , 
Te = ’ { ca (ih/Ty7h,) Je,om (si, u) S€,0m(S1, V)A€,Om 
m=0 


+ (ius /T17hy,) Je,0m (51, U4) S€,0m(S1, ¥)be,0m} e ™’, 


E, = > { Je,on(s1, u)S€,0m(Si, 2)ae,0m} € ”’, 


m=0 


Hy = & { (iki /Trurwh,) Je,om (si, 4) S€,0m(S1, 0)A,0m 


m=0 
— (ih/T17hy) Je,0m(s1, U) S€,0m(S1, ¥)be,0m} e ””, 
H, = Dd {—(ikv /Tiprwhy) Je,0m(S1, 4) S€,0m ($1, 0)40,0m 
m=0 


— (th/T 17h.) Je,0m(S1, U) S€,0m(S1, ¥)b€,0m} e””, 


o 


H, = YX { Je,om(si, u) S€,0m(S1, 0)b€,0m} €. 


m=0 


The fields are represented as a set of even and odd transverse magnetic and 
transverse electric waves whose amplitudes are given by the coefficients 
ae,o, and be,o, respectively. 4, and h, are the metrical coefficients whose values 
are given by h, = h, = cov/(cosh*u—cos’v). The primes above a cylinder 
function denote differentiation with respect to either the variable u or the 
variable v. 
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By employing eq. (8) a similar set of field equations can be written for the 
external region u > uo. 
PROPAGATION MODES ALONG A COAXIAL LINE 


The problem investigated in this paper can be considered as constituting 
an extension of the general problem of wave propagation along an elliptic 
cylinder (Karbowiak 1954). The method of analysis follows essentially the 
excellent treatment given by Stratton (1941, Chap. 9, Sec. 9.19) to the prob- 
lem of a coaxial cable of circular conductors. Consider two confocal elliptic 
cylinders (1) and (3) whose surfaces are denoted by uo and u; respectively. 
A cross section of the coaxial line is illustrated in Fig. 2. The central conductor 


¥e Ly 










GLE YYV: 
ty Yi 


“a 


Fic. 2. Cross section of coaxial line. 


is assumed to be carrying an axial current. It is required to investigate the 
nature of the propagation modes and to determine the propagation constant 
of these modes. 

The set of expressions in eq. (9) implies that an infinite number of modes of 
propagation are possible along a solid conducting cylinder. However, for the 
case of a conductor carrying an axial current, it can be demonstrated that 
only the transverse magnetic mode corresponding to m = 0 possesses a rela- 
tively low attenuation. The asymmetric modes corresponding to m > 0 are 
highly attenuated, consequently they never play a part in the propagation of 
current along a conductor. 

In practice the outer conductor of any coaxial transmission line is of finite 
thickness. However, for the purposes of this analysis it is necessary to assume 
an outer conductor of infinite thickness. This assumption is justified if we 
assume that the skin depth is very much less than the thickness of the outer 
conducting cylinder. In effect, this means that the analysis is restricted to 
frequencies whereby this condition is valid. It will be seen that in order to 
obtain a solution for the propagation constant, it will be necessary to assume 
that the conductors possess infinite conductivity. 

The analysis requires the determination of the radial impedance function 
at each of the boundary surfaces. The impedance Z, can be found by applying 


the definition 
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(10) = =Zfi,. 


Utilizing the results of eq. (9), one obtains the following expression for region 


(1), 


(11) Mz Be T, ‘urwhy Jeo(ss, %) u) 


ik,’ Jeo(si, u) 


To construct a solution appropriate to region (2) which will satisfy the 
internal and external boundary conditions, two independent Mathieu functions 
must be used. Consequently for region (2), uw <u < um, let us write 


(2) __ tke? 
(12) H, = eo 


where Neo(so, u) is an even radial Mathieu function of the second kind, and 
A and B are undetermined coefficients. 
The radial impedance for region (2) becomes 


(13) (2) 7 ei Tes M2ioly uA Jeo(so, u)+BNeo(so, u. u) 
e ike’A Jeo(Se, u)+BNeo(so,u) 


AA Jeo(s2, u) Seo(s2, 0) +BNeo(se, u)Seo(se, a} oo 


In the outer conductor, the proper behavior of the field at infinity is ensured 
by choosing a Mathieu function of the fourth kind. Hence for u > 1, let 


(14) OH, = ft cot (s3, u) Seo(ss, eben 


where C is likewise an undetermined coefficient. 
The radial impedance in region (3) becomes 


—P3"uawhy Fe? (53, u) 
iks Heo (S3,u) 





(15) Zu = 
The coefficients A, B, and C are determined by applying the continuity 
condition of the radial impedance functions, namely: 
ia <- ey. at u = Uo, 


(16) 
ee 


ll 


3. 
OZ, atu=u. 


The substitution of eqs. (11), (13), and (15) in the above relations yields the 
following set of equations: 


rT; Ti'uwwhy Jeo(ss, uo) uo)  T2'uowh,A Jeo(S2, uo) +BNeo(se2, uo) 





tk; 1 Jeo(s1, Uo) - tk: A Jeo(So, uy) + BNeo(se, Uo) 
(17) Ps‘usahylfeo (ss, 1) ao Te uowhy vA Jeo(ss, ui) +BNeo(Se, u1) 
ik;'Heo ( (S3, U1) ike"A Jeo(so, Ui) +BNeo(S2, 1) ~ 


Rearranging the terms in eq. (17) and solving for the ratio A/B, one obtains 
finally 
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QT; ry wake *Jeo(ss 1 uo) 
_T,’ 2 Mek *Jeo(ss, uo) 


Trimk “uke” *Teo(s1, uo) uo) : 
To" uoki *Jeo(ss, Uo) Jeo(ss, we) 


Neo (Se, me? = 


Jeo(so, Uo) — 


Neo(so, Uo) 





a) 8 - 


by] ae 


-3ee *uske "He? (Ss, 1 (Sa 1) 7, 
Neo(ss, u1) ‘ei makes Hey (Ss, (ss, U1) | Neo(ss, 11) 
ie Tusk ‘Hes (S3, U1) 


rT, “wos "Hes (sz, Uy) 


The roots of this determinantal equation are the characteristic values 
ho (pb = 1,2,3,...) which establish the allowed symmetric modes of 
propagation. 

However, it is necessary to resort to an approximate method for determining 
the roots of eq. (18), by assuming that the inner and outer conductors have 
infinite conductivity. From eq. (6), it is recalled that [,? = k,?—h? and 

= k3;?—h?. The imaginary part of kh must remain finite if a wave is to be 
propagated. Hence as the conductivity approaches infinity, T;—>%, and 
I; — k;, and at the same time the absolute values of both k; and k; approach 
infinity. Under these conditions eq. (18) reduces to 
A _ Neo(s2, uo) _ Neo(se, u1) 


C <a oe Nieman ee 
(19) O” Taine Jaw 


Eq. (19) is satisfied when sz = 0. By definition, sz = co?I'2? and since ¢p is 
always finite for all practical cases, it follows that the principal root of eq. 
(19) corresponds to I, = 0, or 


(20) h = ko. 


In addition to the principal wave, higher-order modes may also exist. If Tz 
is large, then the following asymptotic expressions can be employed (Stratton 
1941, Chap. 6, Sec. 6.12), 


Jeo(s2, U1) — Jeo (so, 1) 














1 
V (cole I’, cosh uo) 
1 
V (co rT; cosh 19) Uo) 


Hence eq. (19) simplifies to 


sin(co Tz cosh wo—(4/4)) _ sin(co T2 cosh u;— (4 /4)) 
cos(¢o I, cosh up—(2/4))  cos(co I's cosh uy— (4 /4)) ° 


Jeo(s2, uo) > cos(co T'2 cosh u»— (2/4)), 
(21) 


Neo(s2, uo) > = sin (co 2 cosh u9— (9 /4)). 





(22) 
Simplification of eq. (22) yields 
(23) sin Co T'2(cosh u,;—cosh uo) = 


Recalling that T,? = k2?—h?, eq. (23) yields an expression for the propagation 
constant of the higher-order symmetric modes. There results 


” 9D cc ccigiecllll esac m 
(24) ho = W/ & F= u,—cosh =| Pe BES Bin + od 
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In eq. (24) note that co cosh vu; = ay, the semimajor axis of the outer conductor, 
and co cosh uo = do, the semimajor axis of the inner conductor. It is apparent 
that the condition for the propagation of the higher-order modes is a function 
of the relative sizes of the two conductors. Further, by equating the two 
terms under the radical sign of eq. (24), one obtains a simple expression for the 
critical frequencies of the higher-order modes. 

Consider again the condition for the propagation of the principal wave. 
For [, — 0, and utilizing the equivalent series expansions for the functions 
Jé(s2, u) and Néo(s2, u), it can be shown that in the limiting case Jéo(s2, u) > 0 
and Neo(se, u) > »/(2/mr) tanh u. Consequently in eq. (12), A = 0 and the 
expression for the magnetic field in region (2) becomes 
25 ey ae sim fa s. ~ thy 
(25) A, = 3 [BV (2/m) tanh u Seo(se, v)} e& ’. 
rT. powh, 

The longitudinal current can be obtained from Ampere’s law 
(26) I= £°H,h, dv 
where the current is of the form J = Iye~*?*. Hence, 
T2pow!o 
27 se 
ot) 2/2 .i ko tanh u 
Inserting B in eq. (25) yields 
( 28) OT, = 1/29 hy. 


The transverse voltage along the line can be found by integrating the radial 
component of the electric field from uo to u, 
wpol 


uo 
(29) V= f EE, h, du = oat, (uy— Uo). 
uy 


Finally, the characteristic impedance Zo, which is defined as the ratio of the 
transverse voltage to the longitudinal current, is given by: 





(30) Zo = 52g, (ui uo). 


For the special case of the inner conductor consisting of a flat strip, uw) = 0; 
hence the characteristic impedance becomes simply 


(31) Zo = ae cosh™' (a;/co). 


In passing it should be mentioned that the characteristic impedance can also 
be determined from electrostatic considerations. A solution using a conformal 
transformation method is given as an example in Smythe (1939). 


CONCLUDING REMARKS 


A formula has been derived for the propagation constant of the allowed 
symmetric modes of propagation. It was discovered that a principal wave 
exists whose propagation constant is that of the medium bounded by the two 
conducting cylinders. Higher-order modes may also exist, but in a practical 





WONG: COAXIAL TRANSMISSION LINE obol 


line they are characterized by high attenuation and therefore belong to the 
class of evanescent modes. The preceding analysis may be useful in providing 
an approximate theory for a coaxial line consisting of a strip inner conductor 
surrounded by an outer sheath of rectangular cross section. 
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ON THE EFFECT OF SPEED ON THE KINETIC FRICTION 
OF SOME PLASTIC MATERIALS ON ICE! 


By C. D. NIVEN 


ABSTRACT 


Curves are shown illustrating the effect of speed on the friction of teflon, 
bakelite, nylon, terylene, plexiglas | (perspex), on acetate, polystyrene, and 
polyethylene, sliding on ice at —5° C. and —15° C. Bakelite and teflon had the 
lowest kinetic friction of the materials Soe Terylene also slid easily but 
polyethylene had rather high friction. Departure from Amontons’ Law was very 
evident at high speed. 


INTRODUCTION 


Very little is known about the friction of plastic materials on ice and the 
work described below was originally started with a view to tabulating friction 
constants on ice like the other physical constants of plastics. The investi- 
gation soon revealed that friction on ice depends on several factors besides 
temperature. While the work was in progress a paper was published by 
Bowden (1953) in which it was concluded that teflon (polytetrafluoroethylene) 
had a much lower static friction than perspex (methyl methacrylate), nylon, 
or terylene (ethylene terephthalate) on packed snow, and that Amontons’ Law 
was obeyed for static friction values. 

The writer’s results had, however, shown up to that date that kinetic friction 
on solid ice violated Amontons’ Law and was definitely dependent on loading. 
Bowden gave an order of magnitude for the kinetic friction values of some of 
the plastics but did not show drag vs. load curves for moderate or higher 
speeds. The work described in the present communication concerns the effect 
of speed on the drag values of samples 10 cm.? loaded heavily. 


METHODS 


The work was done on the turntable described already (Niven 1954); the 
method of loading has also been described and accounts for the non-integral 
load values of 14.9, 39.1, 63.3, and 87.6 kgm. being used in the measurements. 
The drag was measured as before by direct weighing. The table was driven by 
means of a d-c. motor, the speed of which could be varied by a resistance. It 
was found that when the speed was cut down the motor had insufficient power 
to drive the table owing to the ‘‘stick-slip’’ phenomenon so well known in 
friction measurements. This trouble of course entirely disappeared at high 
speeds, but at very low speed the table had to be turned by hand or the motor 
helped by hand. The exact speed of revolution was thus not known accurately 
and is simply designated as under 12 r.p.m. Another set of measurements was 
made at about 20 r.p.m. and a third set at around 80 r.p.m. These speeds 
correspond to “less than 1.1 m.p.h.’’, 1.8 m.p.h., and 7.1 m.p.h. However in 
view of the fact that the slider was sliding in the same track at each revolution, 
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the drag values are probably too low--particularly at high speeds—and in 
order to avoid the possibility that the casual reader might take too much out 
of the graphs, r.p.m. values and not m.p.h. values have been assigned to the 
various curves. 

The results must be considered approximate; in most cases at least two 
points are given for drag values, but in many instances more points were 
available for plotting. Getting exact repetition in friction measurements is a 
constant source of disappointment to the experimentalist and when conditions 
are such that the slider starts the stick-slip motion it becomes almost guess- 
work to assign a drag value at all. Poor repetition at slow speeds and low 
temperatures was therefore to be expected. At high speed repetition was, 
however, good enough to permit of a satisfactory estimate being made of the 
position of the drag vs. load curve. To indicate how often the same value for 
a drag measurement was obtained arabic numerals have been attached to 
repeated points, and so in tracing the curves points with a high numeral 
attached have been favored. These curves are approximate and have been 
drawn mainly to show at a glance the important effect of speed. 

The samples consisted of pieces of material 2} cm. wide and 43 cm. long, 
which gave an area of 10 cm.? after the front edge had been rounded off. Since 
perspex, lucite, and plexiglas are trade names for the same material, the 
chemical name, methyl methacrylate, has been used to avoid ambiguity. 
The material called terylene was, according to the supplier, “‘essentially 
polyethylene terephthalate with an admixture of titanium oxide as a delustrant. 
It is made by the addition reaction of ethylene glycol with terephthalic acid 
or an equivalent compound.”’ Pliofilm is a rubber hydrochloride and was 
received from the manufacturer as a thin sheet 0.0025 in. thick. The vinylite 
sample was cut from a flexible sheet 0.036 in. thick. These two samples had 
to be specially clamped to small blocks of smooth flat brass for the tests. 
Since Bowers, Clinton, and Zisman (1954) have found that a thin film of a 
polymer on a metal surface has lower friction than the bulk polymer, it is 
unsafe to compare the curves for these two samples, particularly that for the 
thin pliofilm, with those for the other samples. They are shown in Fig. 5 but 
will not be referred to further in this article. 


RESULTS AND DISCUSSION 


The results are shown in Figs. 1-5. Fig. 1 refers entirely to teflon; so many 
observations were made on this important material that clarity demanded a 
separate diagram for each speed. Figs. 2, 3, and 4 refer to bakelite, terylene, 
methy! methacrylate, polystyrene, cellulose acetate, nylon, and polyethylene. 
They have been grouped roughly according to sliding qualities. Fig. 5 was 
mentioned above. 

The figures show graphically what Bowden mentioned in the text, namely, 
that when sliding speed is increased there is a very marked drop in friction— 
in other words speed is another variable which has to be taken into account 
in evaluating an ice friction constant. The graphs also show that the drag vs. 
load relation curves toward the load axis at high speed more markedly than 
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Fic. 1. Drag vs. load curves for teflon at —5° C. and —15°C. 
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Fic. 2. Drag vs. load curves for bakelite and terylene at —5° C. and —15°C. 


at low speed. Thus there is a suggestion that as speed is lowered to zero the 
curvature may disappear altogether and the relation become linear. Bowden 
working with loadings of less than 6 kgm. on 75.4 cm.? found a linear relation 
for velocities close to zero; such loadings correspond to about 1.2 kgm. on the 
writer’s graphs. Over such a small range departure from the linear relation 
might be quite difficult to detect; nonetheless the points for perspex and nylon 
given in Fig. 7 of Bowden’s paper would permit of a slight curvature for the 
relation which has been drawn in that diagram as a straight line. 

The finding that ice friction decreases with increase of speed fits in quite 
well with the conception that the tips of the asperities melt by ‘frictional 
heat’’. This being so, the usefulness of the turntable in making precise measure- 
ments of friction may be limited because the frictional heating effect must be 
cumulative, when the slider passes over the same track again and again. On 
the other hand the fact that the friction drops immediately as soon as the 
turntable is speeded up shows that the cumulative effect is secondary and not 
the main cause of the reduction in friction with speed. 
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Fic. 4. Drag vs. load curves for nylon and polyethylene at —5° C. and —15° C. 


The graphs for the various plastics do not entirely endorse such an out- 
standing superiority for teflon as Bowden implies nor do they permit of 
grouping nylon, terylene, and methyl methacrylate together. Since the curva- 
ture on the writer’s graphs below 1.5 kgm. is not accurately known, it is hard 
to make a comparison with Bowden’s values. Table I must therefore be con- 
sidered as giving merely a rough estimate of the kinetic friction of teflon, 
nylon, terylene, and methyl methacrylate—at —5° C. and —15° C.—loaded 
with less than 1 kgm. on a 10 cm.* area and moving over the ice very slowly. 

Bowden gives a value for “‘real ski" of 0.045 at —5° C. and 0.05 at —11°C. 
for teflon on packed snow. The writer's estimate is not too far out of line with 
these values. However Bowden's values of 0.34 for perspex and 0.30 for nylon 
on packed snow at — 10° C. do not fit in with the writer’s results. Small changes 
in r.p.m. at low speeds do undoubtedly make large changes in friction so there 
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Fic. 5. Drag vs. load curves for pliofilm on brass and vinylite on brass at —5° C. and 
—15°C. 








TABLE I 
Material —5°C. —15°C. 
Teflon 0.03 0.055 
Nylon 0.045 0.1 
Terylene 0.07 0.075 
Methy! methacrylate 0.07 0.12 


is the possibility that the writer’s lowest speed was still too high to simulate 
static friction. 

Bowden's instrumentation did not permit of precise values being given for 
the kinetic friction but his estimate of 0.02 for a kinetic friction value at 
—10° C. fits in with the writer’s results. There is, however, no possibility of 
using the writer’s graphs to establish the fact that teflon can be assigned a 
friction value, under the experimental conditions used, equal to one third or 
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one quarter of that of nylon, perspex, or terylene; nor can the statement be 
endorsed unconditionally that ‘‘the friction is directly proportional to the 
applied load’’. It is untrue at high loadings and high speeds. 

With the present state of our knowledge of the friction of ice in general, it 
is premature to assign friction values to the various plastic materials because 
these values depend not only on temperature but on loading, speed, and the 
condition of the ice; some believe that even the atmospheric humidity plays 
a part. Ice simply does not conform to the ordinary laws of friction and the 
explanation of its anomalous behavior almost certainly depends on the fact 
that it expands on freeezing. One other point, however, must not be over- 
looked—the observations are being made on this anomalous solid near its 
melting point. 

The graphs submitted above do nonetheless give an idea of the order in 
which the various plastics must be arranged for kinetic friction value. It must 
be stressed, however, that this is merely a suggested order and is not based 
on a statistical treatment of a large number of measurements. 

At —15° C. the suggested order is as follows: 1. bakelite (smallest) ; 2. teflon; 
3. terylene; 4-6. methyl methacrylate, cellulose acetate, polystyrene; 7. nylon; 
8. polyethylene (largest). 

At —5°C. it is somewhat harder to arrange the materials but the graphs 
lead one to suggest the following order: 1, 2. teflon, bakelite; 3-7. terylene, 
methyl methacrylate, polystyrene, cellulose acetate, nylon; 8. polyethylene. 

The measurements at —15° C. simply will not allow nylon to be bracketed 
with terylene; other measurements not reported here endorsed that conclusion. 


CONCLUSIONS 


(1) At high speeds on ice the departure from Amontons’ Law is more pro- 
nounced than at low speeds. 

(2) At —15°C. high speed reduces the friction more markedly than at 
—5°C. 

(3) Teflon and bakelite have lower kinetic friction values on ice than most 
other plastics have. Polyethylene has a comparatively high one. 
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END-FIRE ARRAYS OF MAGNETIC LINE SOURCES MOUNTED 
ON A CONDUCTING HALF-PLANE! 


By R. A. Hurp 


ABSTRACT 


Formulae are developed for the far-field radiation patterns of end-fire arrays 
of magnetic line sources mounted on the surface of a half-plane. Patterns have 
been plotted for a number of interesting cases. It is found that the Hansen- 
Woodvard optimum end-fire condition no longer holds when the array is near 
the edge of the half-plane. The theory seems to describe reasonably well the 
behavior of corrugated surface radiators embedded in a finite ground plane, 
providing the array is removed from the edge by a distance about equal to the 
array length. 


INTRODUCTION 

Although there is a considerable body of literature dealing with various 
aspects of diffraction by a half-plane, the problem of antenna arrays in the 
presence of a half-plane seems to have received little attention. Moullin (1949, 
pp. 196-198) has computed some patterns for arrays radiating broadside to 
the half-plane, while Elliott (1954) has considered approximately the effect 
of a finite ground plane on the radiation pattern of a corrugated surface. 

‘In this paper we derive expressions for the radiation patterns of arrays of 
magnetic line sources (e.g. slots) embedded in a perfectly conducting half-plane. 
The sources are taken to be infinitely long, and are presumed to run parallel 
to the edge of the plane (in the z direction), making the problem two-di- 
mensional. It is assumed that the sources radiate only into the region above 
the plane, and that the phasing is such that the array would normally end-fire 
towards the edge of the plate. The effect of the plate is to cause the direction 
of fire to be tilted away from the line of the array, as indicated in Fig. 1. We 
shall be interested in finding out how the beam tilt, the beam width, and the 
side-lobe level behave as functions of the various parameters of the array. 


EXPRESSION FOR THE RADIATED FIELDS 


Suppose, in the (p, ¢, z) co-ordinate system, we have a plane wave impinging 
on the plate from a direction ¢; Fig. 2. The wave is assumed to have only a 
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Fic. 1. An end-fire array mounted on a half plane. 


1Manuscript received December 27, 1955. 
Contribution from Division of Radio and Electrical Engineering, National Research 
Council, Ottawa, Canada. Issued as N.R.C. No. 3886. 
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xX INCOMING WAVE 


Fic. 2. The co-ordinate system. 


z-component of magnetic field. Then from Baker and Copson (1950, pp. 142- 
144) the field at a point (x, 0) on the screen is 


V2kz- cos ¢/2 


(1) H,(x,0) = 5 exp(ia/4+ikx cos $) f e dx. 


A time dependence of e“‘ has been assumed, and k = w/c. 
By the principle of reciprocity Equation (1) gives the radiation pattern of 


a magnetic line source at (x, 0). 
Consider now an array of N+1 line sources on the plate a distance d apart. 


The amplitudes are assumed equal, but the phase of each is taken to be Bd 


relative to the next. The pattern of the array becomes* 
ix/4 N V[2k(D+nd)] cos ¢/2 


(2) H= = 2 exp[isnd+ik(D-+nd)cos ¢) e dx. 


Here D is the distance from the edge to the first element of the array. Making 
use of the integral 


0 -_ 
(3) f e dx = 1/2\/re™", 
Equation (2) becomes 


(4) H,= L {1+ (1+7)[C(y) -—iS(y)]} exp[i6nd+ik(D+nd) cos ¢] 





where y = 2k(D-+nd)cos*¢/2, and where C(y) and S(y) are the Fresnel 
integrals defined by 


C(y) = Be Jone &,. 


1 ae 
S(y) = Van J sin x ge 
In Equation (4), the (+) sign is to be taken in the “lit” region 0 < ¢ < x, 
the (—) in the shadow region r < ¢ < 2r. 


*The case of arrays of magnetic or electric line sources making some arbitrary angle with the 
half-plane could be treated by an easy extension of the present analysis. 
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BY p- 


Equation (4) is cumbersome for computational purposes, especially when 
N is large. A considerable improvement results if we consider a continuous 
array; that is, we let d ~0, N — o, in such a way that Nd remains constant 
and equal to L, the array length. The summation is replaced by an integration 
which can be performed easily. Equation (4) becomes 











ett k+y)u u=D+L 
(5) H, = oe 
oo u=(8—k)(D+L) 
+ a—k a ree [C(yu)+iS(yu)] er 


where y = 2k cos’p/2. Here again the (+) signs pertain to the region 
0 <¢ < 7m, the (—) signs tom < ¢ € 2r. 

Henceforth our calculations will be based on Equation (5). Whether or not 
a continuous array is a good approximation of a discrete one will depend on 
how large d is. As indicated later, the approximation seems good for d as large 


as \/4 at least. 
It is of interest to examine several limiting forms of (5). If we let 8 > k, we 


obtain an end-fire array with free space phasing, whose pattern is given by 


oe 


” aoa 
(6) H,= se (1+i)1C(rw) —éS(m) lhe é—1)4/ | 
| emp 


If we let D become very large and assume that ¢ is sufficiently far removed 
from m so that yD is large, we have C — 1/2, S — 1/2, and therefore 


i(B—k+y)u | u=D+L 





(7) H, = —2i =e Fac for 0 < @< x, and 
(8) H, = O(yD)“* ior r< ¢ < 2z. 


Equation (7) is just the pattern of an end-fire array embedded in an infinite 


ground plane. 
When ¢ = 7, so that y = 0, Equation (5) simplifies to 


e i(8—k)u | u=D+L 


9) H, = ——-~ 
(9) + * 6—B) leno 
THE RADIATION PATTERNS 


When y and D are such that the approximation of Equation (7) does not 
apply, it is necessary to resort to plotting Equation (5). This has been done 


in Figs. 3-7 for different values of (@—&)L in the range 0 to 22, with a = 
D/(D+L) = 0. Normalized power is here plotted against the abscissa n = 


(2kL)!” sin ¥/2, where ¥ = r—@. 
If we exclude the case (@—k)L = 2x, we ‘observe that, for constant and 


large L, and increasing B—k: 
(a) the main beam tilt decreases, 
(6) the beam width increases, 
(c) the side-lobe level increases. 
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POWER 
NORMALIZED POWER 


NORMALIZED 


NORMALIZED POWER 


Fic. 8. Power pattern for 8 = k, a = 
Fic. 9. Power pattern for 8 = k,a@ 


Fic. 10. Power pattern for 8 = k, a = 
Fic. 11. Power pattern for 8 = k, a 
Fic. 12. Power pattern for 8 = k, a = 


WORMALIZEO POWER 


It thus appears that the directivity decreases continuously as 8—k increases. 
This is in direct contrast to what happens in a normal end-fire array, where 
the gain is a maximum for (@—k)L = a (see Hansen and Woodyard 1938). 

The pattern corresponding to (@—k)L = 2m is an interesting one (Fig. 7). 
Apparently here the ‘‘main”’ lobe has been pushed into the shadow region, as 
evidenced by the low maximum near n = —2; and the ‘‘first side-lobe’’ has 
become the new main beam. 

A second set of patterns is shown in Figs. 8-12. Here we take 8 = k and 
let a range between 1/3 and 1. The abscissa in each case is taken to be 
n = V/(2k(D+L)] sin y/2. For D+L constant and large we note that the half- 
power beam width is roughly constant up to a =3/4, at which point the beam 
width becomes meaningless because of high side-lobes. At a = 1/3 the side- 
lobe level is lowest, being comparable to that of an ordinary end-fire array 
without ground plane; but the beam width is narrower than that of the en- 
hanced end-fire array of length D+L with no ground plane. For a = 1 
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(Fig. 12), we get the expected pattern of a point source on a half-plane. Fig. 12 
also gives the pattern for small y of a finite array at a large distance from the 
edge. 

In Figs. 13-16 we plot the interesting case of (8—k)L = m for various 
a values; (the curve for a = 1 is the same as when 8 = k). The abscissae are 
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Fic. 13. Power pattern for (@—k)L = 7, a = 1/3 
Fic. 14. Power pattern for (@—k)L = 1, a = 1/2 
Fic. 15. Power pattern for (@—k)L = 7, a = 3/4. 
Fic. 16. Power pattern for (@-—k)L = x, a = 4/5 


BEAM TILT 7o 
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Fic. 17. The beam tilt 9 as a function of a. 
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again » = +/(2k(D+L)] sin y/2. Here the beam width and the side-lobe level 
reach minima for a between 1/2 and 3/4. Thus the directivity is maximal here. 
In Fig. 17 we show the variation of the beam tilt 79 with a for the cases 8 = k 
and (8—k)L = =. 

Comparing the patterns for 8 = k with those for (@—k)L = 7m it is seen that 
the directivity of the former is greater for a up to about 1/2. For a = 3/4 the 
(8—k)L = m pattern has become definitely more directive. 

ANALOGY WITH SURFACE WAVE ANTENNAE 

It is to be expected that corrugated surface and dielectric slab radiators will 
perform in a manner analogous to the end-fire arrays here considered. It is of 
interest then to compare previous theoretical and experimental results on 
corrugated surface antennae with the present theory. 

In an interesting paper, Elliott (1954) has calculated the effect of a finite 
ground plane on the radiation pattern of a corrugated surface. We compare in 
Figs. 18 and 19 our results with Elliott’s theoretical and experimental results. 


wees ELLIOTT'S THEORETICAL 
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Fic. 18. Comparison of experimental and theoretical results on the beam tilt. 
Fic. 19. Comparison of experimental and theoretical results on the beam width. 


It is seen that insofar as the beam tilt is concerned, Elliott’s theory agrees 
more nearly with experiment. On the other hand, our values of beam width 
agree quite well with experiment.* 

It appears from Fig. 18 that our theory is not too applicable to corrugated 
surfaces as far as the beam tilt is concerned when a is less than 1/2. For a 
greater than 1/2 agreement is reasonable. We have seen that at a = 1/2 the 
directivity is much the same for 8 =k as for (8—k)L = m. It is therefore sug- 


*Elliott does not give a theoretical value for the beam width, other than a lower bound. 
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gested that the end-fire condition (8—k)L = ma may not be optimal for corru- 
gated surfaces with a near 1/2. 
An experimental study of this point seems desirable. 
COMPARISON WITH DISCRETE ANTENNA 


In Fig. 20 we compare the pattern of a discrete array having d =i\/4, 
D=0,L = 2h, (@—k)L = x with the relevant continuous array. The patterns 


NORMALIZED POWER 











Fic. 20. Comparison of discrete and continuous array. 


show generally good agreement, and appear to justify basing our computations 
on the continuous distribution. 
CONCLUSIONS 

Exact expressions have been derived for the far fields of an end-fire system 
of equally spaced magnetic line sources mounted on a half-plane. The continu- 
ous array appears to be a reasonable approximation to a discrete array for 
source spacings of up to \/4. Three sets of patterns have been computed for 
the continuous array; they are: 

(i) ground plane length D = 0, (@—)L ranging between 0 and 2z, 

(ii) 8 = k, various values of a = D/(D+L), 

(iii) (@—k)L = zx, various values of a. 
The continuous end-fire array seems to represent fairly well a corrugated 
surface with finite ground plane provided a@ > 1/2 approximately. 
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THE RELATIONSHIP BETWEEN THE STATISTICAL 
AND FIELD THEORETICAL TREATMENTS OF 
MULTIPLE MESON PRODUCTION! 


By YASuSHI TAKAHASHI? 


ABSTRACT 
The aim of this note is to point out the relation between Fermi’s statistical 
theory of multiple meson production and the field theoretical treatment. 

Since Heisenberg emphasized the theoretical importance of multiple meson 
production as a phenomenon beyond the applicability of quantum field theory, 
many investigations have been made. Among them, the LOW theory and the 
statistical discussion by Fermi are well known (Fermi 1950; Lewis 1951; 
Milburn 1955). These theories are based on assumptions considered to be 
appropriate for fast and slow collisions, respectively. By fast we mean that the 
collision time is very much shorter than the reaction time. 

As is well known, the LOW theory is valid when the Bloch—Nordsieck 
approximation is good. On the other hand Fermi’s idea is that the energy will 
be concentrated in a small volume during the collision and statistical equi- 
librium will be reached so that the probability that a certain number of 
particles are created is essentially determined by statistical laws. This assump- 
tion would be a good approximation if the collision time were long enough 
compared with the reaction time. 

We shall show in this note the relation between these two theories.* 

According to our previous paper (Umezawa et al. 1952), the total probability 
for the production of 2 mesons by a nucleon is 


4 = r 2 

(1) w= tf. Sao on) AN (ki)... AN(k,) . |v(AP, AE)|", 
(con. ) 

where (con.) signifies that the integration has to be carried out in such a way 

as to satisfy the conservation laws 





| 
2) |AP = yk 
az = ex, 


'1Manuscript received October 28, 1955. 
Contribution from the Division of Pure Physics, National Research Council, Ottawa, 


Canada. Issued as N.R.C. No. 3885. 

2National Research Laboratories Postdoctorate Fellow; now at Department of Physics, 
State University of lowa, Iowa City, Iowa. 

8The importance of this problem was specially emphasized by Professor Taketani a few 


years ago. 
‘It is of importance that (1) is derived under the assumption that the O(k) and v(AP, AE) 


(see below) are commutable with one another, i.e., the production of each meson takes place 
in a statistically independent way. _ notation is as follows: 
1g? /A0(k) eee f 
AN(k) = 2K\K—-k VJ (k?+yu2) = K, p» = the meson mass; 
AO(k) = O'(k) —O’'(k) is the difference of O(k) before and after the nucleon scattering by 
the potential v. O(&) is a quantity describing the nucleon state in the meson cloud. 
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If we work with the occupation number representation, (1) can be rewritten 


as 
Wa = 7 a on {ani} na {atvdba)) lv(AP, AE)|? 
1 1 SAN(k,)\™ | i 
= n! » ! 2! anno} . |v(AP, AE)| 
(con.) ” ; 
(3) = ,y dw, (1, %2,...-); 
“sy 
where 
(4) a 1 fanny . |p(AP, AE)|? 
with 
| n= >on, 
(5) J AP = zz k,n,, 


[AE = Kun, 


The quantity dw,(m, m2,...) is interpreted to be the probability for the 
production of m; mesons with momentum k;, mz mesons with momentum kz, 
and so on. 

It should be noted that our observation of multiple meson production 
phenomena is not that of w, but that of the most probable dw,(m, mz, ... ). 
From this point of view, we shall try to find the most probable distribution 
under the conservation condition.’ This is easily done as follows: For the sake 
of simplicity we will assume AN(k) to be independent of the direction of k. 
In this case, the momentum conservation law may be omitted. The most 
probable distribution will be obtained by maximizing dw, under variation 


of all the arguments m, m2, ..., subject to the energy conservation law; thus, 
(6) 2 toe TI 2 > 1 fatty \s "+ ) aE EK} | wt 

This gives 

(7) ns; = ae gm. 

where 6 will be determined by 

(8) AE = > Kym, = (27)~* | KAN(k) € "dk. 


Then the multiplicity » will be given by 
(9) n= >on, = (24) | AN(k) dk. 


5Since we are interested in the production of any number of mesons, the conservation of the 
meson number is not required. 
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We immediately notice the similarity between our equations (8) and (9) 
and Fermi's statistical theory. Our AN plays the same role as Fermi’s ©. In 
fact, if we take Qas AN, which is independent of k but proportional to E~ (E is 
the energy of the incident nucleon), we get from (9) and (8) 

(10) n = 1.42(w—2)*/4/w'/4, (w= E/M), 


which is in agreement with Fermi’s multiplicity. The. above assumption 
exactly corresponds to Fermi’s choice 
2M _ 4x 12M 
— . = 
Our AN(k) has wider freedom than Fermi’s 2 does. 

We may thus conclude that Fermi’s statistical treatment of meson pro- 
duction is not far from the field theoretical discussion, in spite of their dif- 
ference in appearance, and can be considered as a special case where the meson 
cloud spectrum AN(k) is proportional to E-'. 

A detailed account will appear later on. 
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Note added in pot Since this paper was submitted, almost the same discussion of the 
problem by Dr. Ito (to be published in the Progress of Theoretical Physics) was brought to 
my attention. I would like to thank Dr. Ito for informing me of his results before publication. 


*Fermi’s formula (32) has a minor mistake. It should be 1.34(w—2)3/4/w™/4, 




















GAMMA RADIATION FROM THE PROTON BOMBARDMENT 
OF BORON TEN' 


By G. B. CuHapwick,? T. K. ALEXANDER,’ AND J. B. WARREN* 


ABSTRACT 


The gamma rays resulting from the bombardment of B!° with protons of 
energies from 0.5 to 2.0 Mev. have been observed with a sodium iodide scintil- 
lation counter. Capture radiation, of energy 

E, = 8.81+40.05+10/11 E, Mev., 


showed a broad resonance at Ep = 1135+15 kev. At this energy, the radiation 
had an angular distribution of the form 1+(0.50+0.05)cos*@ and a total cross 
section (3.5+1.0)10~* cm.? Several lower energy radiations were also observed 
and assigned tentatively to cascade transitions in C, 

The cross section for the 430 kev. radiation from the reaction B!°(p, ay)Be? was 
found to be 0.21+0.05 barn at E, = 1.52 Mev. This radiation was found to be 
isotropic. 


INTRODUCTION 


The compound nucleus, C'*, formed in the capture by B!° of protons of 
energy less than 3 Mev. decays mainly by alpha-particle emission. Brown ef al. 
(1951) observed two groups of alpha particles corresponding to the formation 
of the residual nucleus, Be’, in its ground and first excited states. The ground 
state group showed two broad resonances, at E, = 1.1 and 1.5 Mev., while the 
other showed only one resonance at E, = 1.5 Mev. Radiation of energy 
9.47 Mev. from the competing reaction B!°(p, y)C" was first observed by 
Walker (1950), who bombarded a thick target enriched in B'® with 1.16 Mev. 
protons and resolved the radiation with a magnetic pair spectrometer. The 
excitation function of this reaction was subsequently studied by two groups 
of investigators with conflicting results. Krone and Seagondollar (1953) re- 
ported resonances at E, = 0.78, 0.95, and 1.33 Mev., while Day and Huus 
(1954) reported only one at E, = 1.2 Mev. and possibly another at E, = 
2.4 Mev. 

In view of this discrepancy we have investigated the capture radiation from 
B'°+ p with a sodium iodide scintillation counter, and have searched for 
possible cascade transitions. A measurement of the cross section for pro- 
duction of 430 kev. radiation from the reaction B!°(p, ay) Be? was also made 
for comparison with previous measurements. 

EXPERIMENTAL 

Targets of magnetically separated B'°, stated to be 400 ugm./cm.? and 
250 ugm./cm.? thick, deposited on gold and platinum, were kindly supplied 
by the Isotopes Division, A.E.R.E., Harwell. These targets were bombarded 
with resolved proton beams of energies from 0.5 to 2.0 Mev. produced by the 


1Manuscript received November 24, 1955. 
Contribution from the Department of Physics, University of British Columbia, Vancouver, 
Be. 
2Now at Cavendish Laboratory, Cambridge, England. 
3Now at Chalk River Laboratories, Atomic Energy of Canada, Chalk River, Ontario, Canada. 
{On leave of absence at National University, Canberra, Australia. 
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University of British Columbia Van de Graaff generator, and were steam 
heated to reduce carbon contamination. Proton energies were determined by 
measuring the accelerating voltage with a generating voltmeter which was 
calibrated with the well-known resonances of the reaction F!%(p, ay)O". 
The gamma radiation was detected with a sodium iodide scintillation 
counter consisting of a 2 in. long by 12 in. diameter sodium iodide block 
mounted on an R.C.A. 6342 photomultiplier. The pulse height distribution 
from the counter was displayed on a 30 channel pulse height analyzer. Gamma- 
ray energies were determined by calibrating the system with the 2.62 Mev. 
radiation from Th C” and the 6.13 Mev. radiation from the reaction F!*%(p, 
ay)O'’. An occasional further calibration was made with the 9.17 Mev. gamma 
ray from the reaction C!¥(p, y)N'4 at Ep = 1.76 Mev. (Woodbury e¢ al. 1953). 
Preliminary runs showed that large gain shifts occurred in the photomulti- 
plier when the counting rate arising from the 430 kev. radiation from the pro- 
lific reaction B'°(p, ay) Be’ became too large. Such behavior has been reported 
previously by Caldwell and Turner (1954). In order to ensure stability, the 
rate of counting pulses representing over 0.4 Mev. energy release in the 
phosphor was monitored separately and kept below 3000 per second. This 
limitation made the time required to obtain a spectrum of the 9 Mev. radiation 
rather long and necessitated long term stability from the electronic system. 


RESULTS 


(a) The 9 Mev. Radiation 
Figure 1 shows a typical pulse height distribution from the counter due to 
the radiation from the reaction B'°(p, y)C"™ taken at 1.0 Mev. proton energy. 
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Fic. 1. Pulse height distribution from NalI(TI) scintillation counter showing the 9.64 Mev. 
radiation from the reaction B!°(p, y)C" at Ey = 1.0 Mev. Counter 10 cm. from target at 90° 
to incident beam direction. 
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The presence of 9.64+.07 Mev. radiation is obvious; the rise in the region of 
7.9 Mev. is not attributed to a separate gamma ray, but rather to differential 
absorption of low energy bremsstrahlung quanta from fast pair electrons 
produced in the crystal. 

Such spectra were obtained for proton energies from 0.5 to 1.6 Mev. in 
100 kev. steps and with the counter at 90° to the beam. The gamma-ray energy 
varied with proton energy in the manner expected for a simple capture pro- 
cess, i.e. 

E, = Q+10/11 E,, 


where E, was taken to be the incident proton energy less one half that lost in 
the target. The proton energy loss in the target was assumed to be the same 
as that for an equal mass per sq. cm. of beryllium, for which element the 
stopping power has been measured by Madsen (1953). Thus, for example, it 
was assumed that a 400 ugm./cm.? target placed at 45° to the beam would 
reduce the energy of 1 Mev. protons by about 120 kev. 

The mean Q value so determined from 10 measurements was 8.81 Mev. 
with a statistical deviation of 0.02 Mev. However in view of experimental 
uncertainties such as target thickness and multiplier drift we would only 
estimate our Q value error to be within +100 kev. The currently accepted mass 
values and p—n mass difference give Q = 8.697 Mev. (Ajzenberg and Lauritsen 
1955), which is perhaps a little low but hardly outside the experimental un- 
certainties. 

The excitation function of the 9 Mev. radiation is shown in Fig. 2. In 
plotting this curve the proton energy has been corrected for target thickness 
as described before. The yield was measured at 90° to the incident beam and 
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Fic. 2. Excitation function of the 9 Mev. gamma ray from the reaction B!%(p, y)C™. 
Proton energy corrected for target thickness, total cross section corrected for the angular 
distribution (1+0.5 cos%). 
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all counts falling in an energy interval of 2.4 Mev. below the full gamma-ray 
energy were taken as proportional to the 9 Mev. radiation, and it is believed 
this interval contained very few arising from 6 Mev. components in the radi- 
ation. The background subtracted was mostly from cosmic-ray events, since 
the runs were of long duration. A small beam-dependent background was 
noticed from radiation of energy above 10 Mev. which, if due entirely to the 
reaction B'(p, y)C™, showed that there was less than 1 part in 300 of B"™ in 
our targets. An angular distribution of the form (1+0.5 cos? 6) has been as- 
sumed in correcting the 90° vield to total yield. 

Because of the general similarity of this excitation function with that found 
by Brown ef al. (1951) for the ground state alpha particles from B!°(p, «) Be’, 
the broad resonance at Ey, = 1.14 Mev. is believed to arise from the same 
level in C" at an excitation of 9.74 Mev. taking Q for the reaction as 8.70 Mev. 
After subtracting an assumed smoothly rising contribution from non-resonant 
capture or from a higher resonance, as shown by the dotted curve in Fig. 2, 
the full width at half height of the resonance curve is 540 kev. with an estimated 
uncertainty of +40 kev. 

The absolute total cross section, o, at the resonance energy was estimated 
to be (3.5+1.0)10-*° cm.*? which is a factor of about two smaller than that 
measured by Day and Huus (1954). This estimate was made by summing all 
counts corresponding to an energy release of over 7.4 Mev. in the crystal. The 
gamma-ray interaction cross section data of Davisson and Evans (1952) give 
the efficiency of the counter as 50%; by postulating a spectral distribution 
shape similar to that for the 6 Mev. radiation from the F!°(p, ay)O"* reaction 
it was estimated that 40% of the 9 Mev. quanta interactions released over 
7.4 Mev. energy in the crystal. This counter efficiency figure, together with the 
solid angle of 0.107 steradians and the target thickness as quoted by the 
supplier, gives the value of ¢ quoted above, the main uncertainties arising from 
counter efficiency and target thickness, which we believe should not cause an 
error of more than 30% in the cross section. 

The angular distribution of the 9 Mev. radiation was measured at the 
resonance energy with the crystal face 20 cm. from the target in order to keep 
the 430 kev. gamma-ray flux at a low enough level and to avoid the need of 
making solid angle corrections. The gamma-ray flux was monitored with a 
second fixed scintillation counter, which was biased to accept pulses corre- 
sponding to an energy release in this counter of between 7 and 10 Mev. To 
avoid errors due to gain drift the runs at each angle were bracketed by runs at 
90°. The resulting angular distribution fits the function (1+0.5 cos’@) quite 
well as shown in Fig. 3. 

The ratio of the yields at 0° and 90° was carefully checked and found to be 
1.5+0.05. The distribution showed no large asymmetry about 90° and some 
of the data gave close equality of the intensity at 45° and 135°, but a small 
asymmetry about 90° could not be ruled out owing to the limited range of 
angles permitted by the lead shielding used. 


(6) The Intermediate Energy Radiation 
The gamma-ray spectrum from 1.5 to 10 Mev. is shown in Fig. 4, for a proton 
bombarding energy of 1.1 Mev. 
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Fic. 3. Angular distribution of the 9 Mev. radiation from B!(p, y)C" at Ey = 1.2 Mev. 
The different symbols indicate renormalizations to runs at 90°. The full curve is a plot "of 
(1+0.50 cos*@). 
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Fic. 4. Intermediate energy radiations from the proton bombardment of B”, background 
subtracted. 


The presence of radiations around 4 Mev. and 6 Mev. is clearly evident. 
While the spectrum is complex three gamma rays of energies 6.55, 4.80, and 
4.35 Mev. appear to be resolved. These could arise from cascade transitions 
through the levels in C" at 6.35, 4.77, and 4.23 Mev. (Ajzenberg and Lauritsen 
1955). The prominence at 2.5 Mev. would indicate radiations at around this 
energy also. There was evidence too, at lower bombarding energies, for a 
gamma ray of about 2.9 Mev. 
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It would be necessary to use coincidence techniques to be sure that these 
were actually cascade transitions. 


(c) The Low Energy Radiation 

The very large flux of 430 kev. radiation from the reaction B!°(p, a)Be™* 
was all too evident in this work. The cross section was measured and found to 
be 0.21+0.05 barn at E, = 1.53 Mev. This value and the excitation function 
for this radiation, shown in Fig. 5, were in good accord with the measurements 
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Fic. 5. Excitation function of 430 kev. radiation from B!°(p, «)Be™. 
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Fic. 6. Low energy radiations from the proton bombardment of B'®. Curve A shows a 
spectrum taken at E, = 1.1 Mev., target on platinum backing. Curve B, E, = 1.2 Mev., 


target on gold. Curve C, background from platinum. Curve D, 1.28 Mev. gamma ray from Na”*. 
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of Day and Huus (1954). The angular distribution of this 430 kev. radiation 
was found to be isotropic to within 2% at E, = 1.2 Mev., consistent with the 
assignment of a spin 3 to the first excited state of Be’. 

Figure 6 shows the spectrum of radiations with energies from 0.6 to 1.5 Mev. 

At bombarding energies above 1.2 Mev. 720 kev. radiation could be resolved 
from the background presumably from B!°(p, p’y)B"®. With targets on both 
gold and platinum backings it was possible to resolve weak radiation of about 
850 kev. energy, which might arise from a further cascade in C"™ or perhaps 
from some impurity in the target. 


DISCUSSION 


These results essentially substantiate the findings of Day and Huus (1954) 
that over the proton bombarding energy range 0.4 to 2 Mev. there is a single 
broad resonance in the reaction B'°(p, y)C" at E, = 1.14 Mev. If the 540 kev. 
width of this resonance, I’, neglecting the barrier penetration correction, is 
assumed to be due essentially to the alpha-particle width, I,, of the corre- 
sponding level in C", the width for re-emission of protons, Ip, can be estimated 
from the cross section of the reaction B!°(p, a) Be’. Using the one-level Breit- 
Wigner formula for the cross section at resonance, 

Oa = 4aRW TT, /T?, 


where W is the statistical factor. Inserting the value for o, obtained by 
taking the differential cross section found by Brown et al. (1951) of 16 mb./ste- 
radian and assuming an isotropic distribution for the ground state alpha 
particles, 

WI, ~ 40 kev. 


Since the angular distribution of the gamma radiation is anisotropic, the 
total angular momentum, J, of the compound nuclear state > 3. Hence W, 
the statistical weight term, > 7. To estimate I,, it will be permissible to 
assume W = 1. Then 

Tl, ~ T'—-Tl, ~ 500 kev. 
and 
Tr, = Ty oy/¢q ~ 10 ev. 


Weisskopf's formula for the transition probability (Weisskopf 1951) for 
9.8 Mev. radiation gives for T, for EZ; radiation ~470 ev., for M,; ~20 ev., for 
E, ~0.2 ev., and for Mz, ~0.001 ev., which might suggest that the radiation 
is M;. 

It is difficult to draw any definite conclusions as to the spin and parity of 
the excited state of C" corresponding to the observed wide resonance, because 
of the lack of information regarding the background in Fig. 2. The yield above 
E, = 1.6 Mev. appears to be either from a higher resonance, as suggested by 
Day and Huus at 2.4 Mev., or from direct capture. Our data are inadequate 
to rule out asymmetry about 90° and hence interference. 

The ground state of C" is probably a (}—) state both on the basis of the 
assignment to the mirror nucleus B" and on shell model arguments. If, then, 
the assumption is made that the observed angular distribution is due to a single 
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state of well-defined angular momentum and if possible interference between 
various multipolarity gamma rays is neglected, some assignments for the 
excited state in C"™ can be rejected. Thus values of J > (], + or —) can 
probably be discarded on the basis of either the observed radiation width or 
the failure to observe cos‘@ terms in the distribution. J = (4, + or —) and 
J = (§+ or] +), formed by s-wave protons, would lead to isotropy of the 
gamma radiation. J = (§, + or —) can probably be eliminated since these 
assignments give a distribution of the form 1 + A cos*@ with A < 0.28 which 
is outside our error of measurement. The possibility of J = ($+) formed by 
1 = 2 protons yields A < 0.41 and, moreover, there is the inherent improba- 
bility of d-wave formation when s-wave formation would be possible. The only 
possibility which has not been eliminated is thus J = (§—), formed by p-wave 
protons and decaying by J/, radiation, which agrees satisfactorily with the 
observed pattern, if a contribution of about 10% from the channel spin 3 
is taken.* 

It must be stressed that this conclusion is based on the drastic assumption 
of no interference between states of different J values and no radiation multi- 
pole interference. An awkward objection to this assignment to the C" state 
involved is the apparent absence of a resonance in the yield of short-range 
alpha particles to the first excited state of Be’ at this energy. Such alphas 
would, of course, be inhibited by their low available energy. 
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AN INVESTIGATION OF THE SODIUM-POTASSIUM 
EQUILIBRIUM DIAGRAM! 


By D. K. C. MacDonaLp, W. B. PEARSON, AND Lois T. TOWLE 


ABSTRACT 


The study of fundamental properties such as specific heat and electrical 
resistivity in pure metals has made it very important to know the influence of 
small amounts of impurity. This has led us to undertake a detailed investigation 
of the sodium-potassium equilibrium diagram. Previous work had indicated zero 
solid solubility at each limit (100% Na; 100% K) but the present investigation 
shows in particular that solubility extends to a few per cent from either limit. 


INTRODUCTION 


The measurement of such parameters as specific heat, electrical resistivity, 
can yield a great deal of information about the state of the lattice in a metal. 
From continuous observations over a wide temperature range on pure alkali 
metals in particular, we have found (MacDonald 1953, 1955; Dauphinee 
et al. 1955; cf. also Carpenter 1953) that both the specific heat and electrical 
resistance exhibit an anomalously rapid rise over a temperature range ex- 
tending some 50° to 100° C. below the melting point. The metals used by us 
were very pure by physical standards (such as the residual resistivity as 
T — 0), and this rise has been interpreted as due to significant concentrations 
of thermal defects in the lattice as the melting temperature is approached. 
On this assumption, both parameters permit us to deduce the energy required 
to create such a defect in the lattice (and these estimates should of course 
agree within experimental error); while from the calorimetric data specifically 
the concentration of defects at any temperature can be determined, and from 
the resistive data we can then deduce the scattering cross-section for electrons. 
Knowledge of the energy required to create the defect and of the scattering 
cross-section may then enable one to identify with some certainty the specific 
type of defect (e.g. a vacancy or interstitial atom—cf. e.g. Machlup (1956)). 
Information about the actual defect concentration may be of value in under- 
standing the melting mechanism, since it is thought that melting may be 
“triggered” by some critical concentration of lattice defects. 

Consequently it is of some importance to be reasonably certain that we are 
justified in identifying the anomalous rise in specific heat and resistance in 
this way. Another possible source of such a rise might be the anharmonicity 
of lattice vibrations which increases as the temperature rises (cf. Born 1933) 
and this led Dugdale and MacDonald (1954) to study quantitatively in a 
specific model (a linear chain) the influence of anharmonicity on the thermo- 
dynamic parameters. This investigation led us to believe that the rise could 
not in fact be assigned to this cause. 

A further possibility was that the presence of even a rather small concen- 
tration of an impurity, insoluble in the solid state but soluble in the liquid, 
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might produce such an anomalous rise as the melting point of the ideally pure 
metal is approached, although the predicted temperature-dependence of 
specific heat due to this is not of the right form. Nevertheless, in our sodium 
and potassium,* it was presumed that the significant residual impurity in 
each case would be the other metal, and since previous determinations of the 
potassium-sodium equilibrium diagram (van Rossen Hoogendijk 1912; Janecke 
1928; Rinck 1933) had all indicated that there was no terminal solid solubility 
at either limit (100% Na; 100% K), it was felt necessary to re-examine the 
equilibrium diagram carefully. In fact the published zero solubility of potas- 
sium in sodium was already contradicted by the observation (cf. MacDonald 
and Pearson 1954) that the residual resistance of initially very pure sodium 
(say Raox./ Re730K, 2X 10-*) could readily be increased by a factor of several 
hundred by alloying with a small percentage of potassium indicating very 
strongly that a significant fraction of the potassium had entered into (random) 
solid solution. 


EXPERIMENTAL METHOD 


A rather extended series of measurements of the electrical resistance of 
potassium-sodium alloys was undertaken and has enabled us to construct an 
approximate equilibrium diagram for this binary system in considerably more 
detail than has previously been available. We may mention that a number of 
our early measurements were undertaken for purposes other than that of 
determining the equilibrium diagram. (e.g. an interest in the resistance of 
liquid alloys for which as yet no satisfactory theory exists), and the importance 
of making measurements on alloys essentially in thermal equilibrium—particu- 
larly in the vicinity of phase boundaries—and also the influence of supercooling 
became clear as the investigation proceeded. The method used of measuring 
resistance as a function of temperature offers considerable advantages in 
investigating an alkali metal alloy system. This is because of the small amount 
of alloy required and also because it is only necessary to make measurements 
on alloys which are strictly in equilibrium when in the vicinity of phase 
boundaries, and these boundaries can be located if necessary in a rapid pre- 
liminary survey. By contrast thermal analysist involves problems of estab- 
lishing steady heating or cooling rates and does not generally allow measure- 
ments to be made on alloys annealed exactly to equilibrium at a particular 
temperature. 

Alloys were prepared from sodium and potassium of rather high purity, 
the residual resistance ratio of the primary metals being 210-4 and ~20 10-4 
respectively. The alloys were prepared in vacuo in a glass apparatus (sketched 
in Fig. 1) which was heated in an oil bath. The molten metals were mixed in 
the U-tube, A, by vigorous shaking and then forced into the mold, B, and 





*These metals were exceedingly pure. For the sodium we are indebted to Messrs. Philips, 
Eindhoven, Netherlands, and Mitcham, Surrey, England, while the potassium was specially 
prepared by I.C.I., Ltd., for the Pure Metals Research Committee of the United Kingdom. 

tTime-temperature curves taken under steady conditions of heating or cooling of the sur- 
rounding heat bath. See, for instance, ‘Metallurgical Equilibrium Diagrams” (Hume-Rothery, 
Christian, and Pearson 1952). 
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Fic. 1. Schematic drawing of glass apparatus for preparing resistance specimens. 


reservoir tube, C, by applying a pressure of one atmosphere of helium after 
the apparatus was inverted. After cooling, the mold was cut free at X and 
sealed with Apiezon compound Q covered with Glyptal, which made an air- 
tight seal that did not exert any appreciable pressure on the metal as it ex- 
panded on heating and melting. 

The filled molds were placed in a simple thermostat and resistance measure- 
ments were made successively with increasing or decreasing temperatures using 
a galvanometer amplifier (MacDonald 1947). Measurements were made gen- 
erally only when the alloy attained equilibrium at the holding temperature, 
particular care being taken about the peritectic and eutectic temperatures in 
the vicinity of the solid solubility limits. In the cooling cycles, however, 
supercooling was frequently encountered at the liquidus and isothermal 
transitions. As no process of seeding was possible the phase boundaries were 
determined only from measurements made during the heating cycles. 

The composition of all alloys examined was determined by chemical analysis. 
When possible both Na and K were separately determined. In dilute alloys, 
spectroscopic analysis was also used to obtain the concentration of the minor 
component. The scatter to be noted in the liquidus points of several alloys 
may be attributed either to the difficulties of sufficiently accurate chemical 
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analysis or to inhomogeneity along the length of the specimens, despite the 
vigorous mixing given to the molten metals. Such inhomogeneity may result 
from the necessity of slowly cooling the molten alloys during preparation in 
order to maintain a continuous thread of metal in the central capillary of 
the specimen molds. Thus one alloy in particular which had an analyzed 
composition of 93.5 at. % K appears from the liquidus temperature and size 
of eutectic arrest to contain effectively about 90 at. % K. 


EXPERIMENTAL RESULTS 


The liquidus curves in the equilibrium diagram shown in Fig. 2 agree rather 
exactly with those obtained by Rinck (1933) by thermal analysis, excepting 


TEMPERATURE (°K) 
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Fic. 2. Sodium-potassium equilibrium diagram. 
X Phase boundaries obtained from resistance-temperature measurements. 
+ Eutectic reaction observed in alloys not completely in equilibrium. 
Quenching experiments in the sodium solid solution indicate that the solid solubility of po- 
tassium decreases considerably at temperatures below 280° C. as sketched. No data are yet 
available for the potassium solid solution. 


that the potassium-rich hypoeutectic liquidus curve lies a few degrees lower 
in the range 80-90 at. % K according to our measurements. On the other hand, 
we, for the first time, have specifically looked for and found an appreciable 
range of solid solution formed by both Na and K. 

The peritectic temperature is determined as 280.5+1° K. where liquid 
containing ~-40 at. % K reacts with the sodium solid solution containing ~3 
at. % K to form-Na2K. At the eutectic temperature of 261+1° K., the com- 
pound Na:K and the potassium solid solution containing ~4.5 at. % Na are 
formed from the liquid containing —66.6 at. % K at the eutectic point. The 
larger solid solubility of Na in K than of K in Na is consistent with the expec- 
tation that it will be easier for a smaller atom to enter solid solution in a 














np eae 


ian ait is 


eS ee ey Ee eee 


MacDONALD ET AL.: EQUILIBRIUM DIAGRAM 305 


solvent whose atoms have a larger diameter than vice versa, particularly when 
the radius ratios differ by 24% as in the present case. 

The terminal solid solubility limits at the peritectic and eutectic tempera- 
tures were determined by bracketing with alloys which did,—and did not—, 
show a discontinuity of resistance as a function of temperature at the invariant 
point. Thus in Na-rich alloys the alloy at 3.1 at. “% K showed no jump in 
resistance while those at 3.7 and 4.5 at. °% K gave a small but definite dis- 
continuous increase in resistance at 280° K. on heating. It was found important 
to ensure that equilibrium was attained in the alloy at each measuring tempera- 
ture in the vicinity of the isothermal transition as incorrect results could 
otherwise be obtained apparently indicating a smaller solid solubility than was 
actually the case. The position of the solid solubility limit was also checked 
approximately by examining the size of the resistance discontinuity, suitably 
normalized, as a function of the alloy composition. Owing to the steepness of 
the solidus curves, bordering the Na and K solid solutions, no accurate estima- 
tion of the solidus temperature was attempted except in the most dilute alloys. 

The composition of the intermediate phase, proposed as Na2K in earlier work, 
was confirmed subsequently by the crystal structure determination of Laves 
and Wailbaum (1942), who showed that it had a hexagonal Cl4 type of 
structure. The structure, therefore, belongs to the group of Laves phases which 
have a composition AsB and owe their stability to the ease with which space 
can be filled with atoms whose radius ratio 7g/r, is 1.225, which is indeed 
close to the value of rx /rx, = 1.24 as determined from lattice-spacing measure- 
ments. The indication that the eutectic horizontal extended to alloys with 
lower potassium content than 33.3% (i.e. beyond the compound NazK) as 
found by van Rossen Hoogendijk (1912), by Rinck (1933), and by ourselves 
is due to lack of equilibrium and incomplete peritectic reaction. van Rossen 
Hoogendijk showed that the thermal arrests at 261° K. in this composition 
range below 33% K disappeared after suitable annealing of the specimens. In 
the present work a discontinuity of resistance could still be found at 261° K. 
in alloys of 26%, 28.5%, and 32.8% (but not in alloy 24.7% (atom K)) even 
after they were annealed at 3° to 5° C. below the peritectic temperature for 
more than eight hours. However, these discontinuities were small and showed 
no regular variation with the composition of the alloy. We are satisfied that 
they were the result of incomplete peritectic reaction during the annealing 
treatment given to the alloys. The present investigation does not therefore 
dispute the accepted composition of the intermediate phase. 
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ELECTRICAL BREAKDOWN IN ARGON AT 
ULTRAHIGH FREQUENCIES' 


By A. D. MAcDONALD ANbD J. H. MATTHEWs? 


ABSTRACT 


Measured values of breakdown electric fields in pure argon gas are presented. 
Measurements were made in two resonant cavities at a frequency of 2800 Mc./sec. 
and for pressures varying from 4X10? to 200 mm. of mercury. The present 
results are in agreement with those of Krasik, Alpert, and McCoubrey. 


In recent years breakdown electric fields have been measured for a number 
of gases at frequencies in the neighborhood of 3000 Mc./sec. (Herlin and Brown 
1948; MacDonald and Brown 1949a, 6). These microwave measurements are 
particularly useful in giving us information about the fundamental collision 
processes in the gas. There are two reasons for this. First, secondary effects 
at the walls of the container are small. Second, with reasonable gas pressure 
variations the phenomena vary from those in which the electron collision 
frequency is much greater than the field frequency to those in which it is 
much less. 

This paper presents measured values of breakdown fields for pure argon. 
Krasik, Alpert, and McCoubrey (1949) made some measurements of break- 
down and maintaining fields at 2950 Mc./sec. in argon, but the present work 
covers a wider range of variation of the experimental parameters. The experi- 
mental procedures are similar to those described previously (MacDonald and 
Brown 1949a, b; MacDonald and Betts 1952). A block diagram of the appa- 
ratus is shown in Fig. 1. Power generated by a c-w. magnetron operating at 
2800Mc. /sec. is coupled to a resonant cavity by a coaxial transmission line. 









CAVITY 
WAVE METER 
BOLOMETER 
MATCHED (THERMISTOR) 


PAD 
| 
GM. 
MAGNETRON 
SECTION on 


Fic. 1. Block diagram of the microwave apparatus. 


A known fraction of the power is measured by a thermistor. The power ab- 
sorbed by the cavity, the cavity Q, and the known field configuration are 
used in the usual manner to determine the electric field (Brown and Rose 1952; 
Rose and Brown 1952). The cavities in which breakdown takes place are made 
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of oxygen-free high-conductivity copper and connected through kovar to an 
all-glass vacuum system. The vacuum system maintains a pressure of less 
than 5X10-* mm. of mercury with the pumps shut off for a period sufficiently 
long to make a series of measurements. All argon used had an impurity content 
of less than one part in 104. Mass spectroscopic analysis by the manufacturer 
indicated that several of the samples used had considerably less than one part 
in 10° impurity. The amount of the impurities present had no measurable 
effect on the breakdown fields. The impurities present were nitrogen, some 
easily removable hydrocarbons, and traces of hydrogen. Ionizations resulting 
from collisions of atoms in low-lying metastable states with impurity atoms, 
which greatly affect breakdown fields in helium and neon, are not important. 
The reason is that the metastable levels in argon are at approximately 11.5 volts 
and are therefore lower than the ionization potentials of the impurities present. 

The measurements were made in cylindrical cavities having heights of 
0.159 cm. and 0.476 cm., and characteristic diffusion lengths of 0.0505 cm. 
and 0.151 cm. The experimental data are presented in Fig. 2, which also in- 


3 
9 
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Electric field (Volts/cm) 





O.! 1.0 10 100 
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Fic. 2. Experimental values of the breakdown electric field at 2800 Mc./sec. for cavities 
having characteristic diffusion lengths of 0.0505 cm. and 0.151 cm. The dashed line represents 
our calculations of electric fields from the data of Krasik, Alpert, and McCoubrey. 


cludes breakdown data calculated from the Krasik, Alpert, and McCoubrey 
(1949) data. The characteristic diffusion length, A, of the cavity they used was 
0.135 cm. The small difference between our data for A = 0.151 cm. and their 
data is what one would expect, and our work is therefore in good agreement 
with theirs. It will be noted in Fig. 2 that at pressures of the order of 0.1 mm. 
of mercury there is a change in the slope of the curves. At these pressures the 
electron mean free path becomes comparable with the cavity dimensions, 
effects at the container walls become important, and the discharge is no longer 
diffusion controlled as it is at the higher pressures. 

The curves in Fig. 2 represent the average of several sets of data and the 
maximum error in the electric field is less than 5% and in the pressures less 
than 2%. 

In previous studies of breakdown in helium, hydrogen, and neon (Mac- 
Donald and Brown 1949a, 6; MacDonald and Betts 1952), the measured break- 
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down fields have been compared with those predicted on the basis of kinetic 
theory calculations using distribution functions derived from the Boltzmann 
equation. Such a derivation requires the solution of a differential equation, 
which includes among its parameters an analytical expression for the cross 
section for collisions between electrons and atoms of the gas under consider- 
ation. For argon this cross section is a very complex function of electron energy 
and a suitable approximation which permits solution of the differential 
equation has not yet been made. 
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REFLECTION AND TRANSMISSION AT A SLOTTED 
DIELECTRIC INTERFACE! 


By R. E. CoLuin 


ABSTRACT 


A theoretical analysis of a slotted dielectric interface is given for the czse of a 
plane wave incident at an angle 6; and polarized with the electric vector parallel 
to the slots. A variational solution yielding upper and lower bounds to the 
parameters of the equivalent circuit for the interface is given. A numerical 
example for a wave incident at an angle of 30° on a slotted dielectric with a 
relative dielectric constant of 2.56 at a wavelength of 3.14 cm. shows that the 
phase shifts introduced into the reflected and transmitted waves are of the order 
of 1°. Also the impedance transformation across the interface is just that due to 
the change in wave impedance. Such a slotted dielectric may be| used to match a 
microwave lens to free space. A theoretical calculation shows that the reflection 
coefficient does not exceed 5% for angles of incidence up to 30° in the wavelength 
range 3-3.28 cm. for a slotted section used to match a dielectric of relative 
dielectric constant 2.56 to free space. Without a matching section the reflection 
coefficient is 23% or more. 


1. INTRODUCTION 


When a uniform T.E.M. wave is incident at an angle 6; on a plane dielectric 
sheet a fraction of the incident radiation is reflected and the rest is transmitted 
into ‘the dielectric sheet and propagates at an angle 6,. For a perpendicular 
polarized wave the reflection coefficient is 


sec 0,— x sec 4; 


* sec 6, + V« sec 0; 


where x is the relative dielectric constant of the material. The reflected wave 
can be eliminated by means of an intermediate dielectric sheet with an effective 
thickness of one-quarter wavelength and a wave impedance equal to the 
geometric mean of those of the two media to be matched. The difficulty in 
producing such a matching layer is due to the lack of suitable materials with 
the required dielectric constants. An equivalent matching layer can be made 
by slotting the dielectric interface. The thickness of the slots is chosen to obtain 
the required wave impedance and the depth is determined so that the section 
is effectively a quarter wave long. In order to design such a matching section 
it is necessary to know the parameters of the equivalent transmission line 
circuit for the interface. One particular form of slotted interface with a perpen- 
dicular polarized wave incident at an arbitrary angle 6; will be considered 
here. The slots are assumed parallel to the incident electric vector. The par- 
ticular case of normal incidence has been analyzed before (Collin 1953; Collin 
and Brown 19550), and the main purpose of this article is to extend the analysis 
to the case of oblique incidence. 

Such slotted surfaces are of considerable value in matching the surfaces of 
dielectric microwave lenses to free space. 
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2. FORMULATION OF PROBLEM 
The type of slotted interface to be analyzed is illustrated in Fig. 1. Each 
dielectric tongue is of thickness ¢ and is separated by an air space of thickness a. 
The spacing of the slots is s and is uniform for all x. A T.E.M. wave with the 
electric vector along the y-axis is incident at an angle 6; on the interface from 


> 





Fic. 1. Slotted dielectric matching section. 


the region z < 0. The incident wave has components E,, H,, and H, and no 
other field components are excited by the structure since it is uniform along 
the y direction. Consequently the total field can be derived from the one scalar 
component E, by means of the following equations: 


(1) tt 2, diel we ae 

The scalar function EF, satisfies the following wave equation: 

(2) OE 5 9 Ee ae) ko Ey = 0, 

where Ro? = wpe = 42r?/Xo? 

and x(x) = « for |x| < }tans,n =0,1,2,...;0<¢2<¢D, 
=<« for 2a 2D, 


= 1 otherwise; 
m.k.s. units are used and the time factor e”' is suppressed. The free space 
wavelength is \o. The boundary conditions to be imposed on Ey, are conti- 
nuity of £, and its normal derivative at all air—-dielectric interfaces. 
The incident electric field will be taken as 
st exp[—jko(x sin 6,+2 cos 6;)]. 


When this wave strikes the interface a reflected wave 
Rs-} exp[—jko(x sin 0,—z cos 6,)] 

plus higher order diffraction modes will be excited in the region z < 0. The 
amplitude of the reflected wave is | R| where R is the complex reflection co- 
efficient. It will be convenient to consider this uniform incident T.E.M. wave 
propagating at an angle @, to the interface normal as a non-uniform wave 
propagating along the z-axis since this leads to a simple transmission line 
equivalent circuit for the structure. 
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In the slotted region waves propagating or attenuated along the z-axis will 
be excited so that F, will be in the form of a summation over normal modes, i.e. 


>> Fm (x) (b4 e77"?-+;, €”"”). 


Substituting into the wave equation yields the following equation for the 
normal mode function F,,: 


8) TEP + (rm +6(x)ho") Fm = 0. 


Because of the periodicity of x(x) a solution of the formT 
Fp = e~* bm (x) 


exists such that ¢, is a periodic function of x of period s (Morse and Feshbach 
1953). Since the incident field varies with x as 
exp (—jkox sin 0;) 


it is necessary for all field components to have this same progressive change 
in phase with x and thus 
(4) B = kosin 0;. 


3. SOLUTION FOR F,,(x) 


Let the two independent solutions to the wave equation (3) in —3s < x < 3s 
be ¥m° and ~»,° having even and odd symmetry about x = 0 respectively. 
(Superscripts e and 9 are used to distinguish between the functions having even 
and odd symmetry.) Since ¥,,°° are solutions to the wave equation, both 
Vn’? and dpm*’/dx are continuous at x = +3. The functions y,,” will have an 
eigenvalue /,, in the region 

|x| < Hans, 2 = 0,1,2,..., 

and k,, in the free space regions. The eigenvalues /, and km satisfy the con- 
ditions 
(5) Yn? = In? — ako? = Rm? —Ro?. 
The solution for all x is F,, and will be a linear combination of ¥,° and Wm’, thus 
(6) Fn = Am¥mn't+Br¥n’- 
For propagation along the x-axis it is necessary that 

Fn(4s) = €7*?*F_(— 4s), 
" Fn! (4s) = "RF q!(— 45), 
where the prime indicates d/dx. Substituting for F,, and F,,’ from Eq. (6) and 


making use of the symmetry properties of ¥,°? allows one to write the fol- 
lowing two equations which must hold for the particular value x = 3s: 


An(1—e7*") Wn? + Bn (1 +e") ym” = 0, 


(8) 
An(1+e-**) pn" +Bm(1—e-**) Ym’? = 0. 


{This is a direct application of Floquet’s theorem. The problem here is similar to the 
quantum-mechanical one of flow of electrons in a periodic square wave potential field. 
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For a solution to exist the determinant must vanish. This condition yields the 
following eigenvalue equation: 
. we bn” 

(9) tan’a(}s) = —Ye Yn) 

Vn Vm’ lre}s 
The eigenvalues /,,, Rm, and ym are determined by this equation when 8 is given. 
However this leads to complicated algebra and a simpler method of deriving 
the eigenvalues by means of matrices is given in Appendix I. 

The coefficient B,, is given by 


(10) By = —Apj tan B(}s) le oe 
The coefficient A,, will be chosen so that 

4s 
(11) Pula’ = 1, 

—}s 


where the star denotes the complex conjugate value. Expressions for A, are 
given in Appendix II. It is readily shown that the functions F,, are orthogonal 
to F,* with m # n. The functions ¢, are given by 

(12) on = 0? F,,. 

Now ¢m(x+s) = e#@+® Fii(x+s) = ec? F(x) = om(x) from Eq. (7). Thus 
om is periodic as it should be. Since ¢n¢,* = F,F,* the functions ¢, are orthog- 


onal to ¢, for m # n also. 
It is now necessary to derive explicit forms for ~,°?. In —}t <x < 31, let 


V.° = cos lis, Vn" = sin Lx. 
In 3t < x < }s, let 
Un = dm? COS km(xX—Om°), Vn? = Ay? Sin Ry (x—8,,°) . 
Then in —3s <x < —it, 
Vimo = Am* COS Rm(X+Om*), = Wn? = Am SiN Rm(X+4Om°) « 
The phase angles 6,,°° and coefficients a,,°° are as yet undetermined. At 
x = +t, yn”? and y,,/“° must be continuous, thus 
COS 1, (3t) = Am* COS Rm(3t—On°), 


Im SiN Lm (3t) = RmOm*? sin ky (3t—-On°), 


and similarly for Ym? and ¥,’°. The coefficient a,° is given by 


(13) Am?’ = COS lm (Zt) SEC Rm (Ft —Om,*) 

and the phase angle @,,° is given by the transcendental equation 
(14) lm tan I_(4t) = Rm tan Rm (3t—On*) « 
Similarly it is found that 

(15) Am? = sin 1_($t) CSC Rm($t—Om’), 

(16) lm COt 1m (Zt) = Rm COt Rm(ZL—Op_°) « 


Eqs. (13) to (16) give solutions for a,,°° and 6,,°° when J, and km are known. 
The equation derived for /,, and k», in Appendix I is 








402 CANADIAN JOURNAL OF PHYSICS. VOL. 34, 1956 






(17) cos Bs = cos [pt cos Rma — fn tin. sin /,¢ sin ka. 
onmm 
Eq. (17) together with the condition ky? = ],?—(k—1)ko? derivable from 
Eq. (5) yields the allowed values of /, and km for given values of ¢, a, and 8. 
The functions F,, are now completely specified and the electric field in the 


slotted section may be taken as 


(18) Do Fm (bme"* +bme™) = 6” 2) bm (bme + bme™*). 
m=1 m=) 











The magnetic field components follow by differentiation as indicated by 
Eqs. (1). The coefficients 5,, will be determined by the continuity conditions 
on E, and H; at z = 0. 


4. SOLUTION FOR THE FIELDS IN REGION z<0 


In the region z < 0 there is an incident wave 
s—} exp[—jko(x sin 6,;+2 cos 6,)] 










and a reflected wave 
Rs~} exp[—jko(x sin 6;—2 cos 6;)] 






plus higher order modes excited. The higher order diffraction field will be a 
superposition of plane waves given by (Berz 1951) 






(19) > 'c, s texp[—jko(x sin 0,—z cos 6,)], 






where the prime indicates the exclusion of the term m = 0. Because of the 
periodicity of the structure and the fact that the field must show a progressive 
phase change along the x-axis it is necessary that 






ko sin 6, = ko sin 0+ 22 , n=+1,+2,..., 







or 
(20) sin 6, = sin ate, nm=+1,+2,.... 
The corresponding value of cos @, is 

(21) cos 6, = /(1—sin’6,) = —jTn ko’. 






When sin’@, > 1 the sign of I, is chosen to give decaying waves along the 
negative z-axis. In order that none of these higher order waves shall propagate 
in the negative z direction it is necessary for sin*@, to be greater than unity 
for all n = 0. From Eq. (20) this imposes the following condition on the 
maximum spacing that can be used: 









Xo 
22 s< 
a) 1+ {sin 6;| 
for all values of Ao and 6; of interest. When such a slotted interface is used as 
a matching section, A» in Eq. (21) must be replaced by Ao/+/x to ensure that 


none of the higher order waves will propagate in the solid dielectric region. 
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This condition also ensures that only one mode will propagate in the slotted 
section, i.e. all y, are real except y:. The total electric field in the region z < 0 
is therefore 


(23) stem | exp (— jes cos 6;) +R exp (jkocos 6;) 


+ 'cnexp|—j(2n 2's ste | ; 
At z = 0, E, and H, or equivalently EZ, and aE, / as must be continuous. 
Thus from Eqs. (18) and (23) 
(24) 14+R+ >! crexp[—j(2nw/s)x] =+/s : (bn +bn) bm 


m= 


(25) (1—R)To- z Cal',exp| —j(2n 4 / dx] = Vs a bn )Y¥mOm 


where Ip = jko cos 6;. The solution to these two equations determines the 
amplitude and phase of the reflected and transmitted waves. A rigorous solution 
would be difficult to obtain so a variational procedure that will give upper 
and lower bounds on the equivalent circuit parameters for the interface will 
be used instead. 

5. VARIATIONAL SOLUTION 


The variational procedure adopted here is a modification of that introduced 
by Schwinger (unpublished). It is fully described by Collin and Brown (1955a) 
in a paper on the junction effect between an empty and inhomogeneously filled 
waveguide. The method is the analytic equivalent of the Weissfloch (1942) 
experimental technique for finding the equivalent circuit of a waveguide dis- 
continuity. The slotted dielectric interface is represented by two transmission 
lines of electrical length u and v connected by an ideal transformer of turns 
ratio n: 1 as illustrated in Fig. 2. The characteristic impedances of the lines 


EZ 
EZ 


Ze= Neo SEC 0; 


ec 


Fic. 2. wae circuit for slotted dielectric interface. 


will be taken equal to the wavelength along the z-axis in the free space region 
and slotted section. With reference to Fig. 2 a field minimum is located a 
distance d in front of the interface when a short-circuiting plane is placed a 
distance / from the interface in the slotted section where 


Zy tan(Bod+u) = —n?Z tan(@,/+02) . 
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The parameters u, v, and » completely specify the effect of the interface on 
the propagating mode and are determined by an analysis of the curve Bod+ {yl 
vs. Bil where 

Bo = Ro cos 8;, b=|n1, 


1 = position of short-circuit from the interface in the slotted section, d = posi- 
tion of an electric field minimum in the region z < 0, measured from the inter- 
face. A unique value for the parameters u, v, and is only obtained provided / 
and d are chosen large enough so that all the evanescent modes have decayed 
to a negligible value at the short-circuit and field minimum positions. 

The variational expressions giving Bod as a function of 6,/ are derived in a 
similar manner as for the waveguide case referred to. For the present problem 


these expressions are 
$s hs i 
—* By cot Bod G dx G"dx 
—ts —hs 


(26) = ‘> 5 c ” Gexpli@n/sirkix fi G*exp|—j(2n2/s)x]dx 


—_ <n 


+, cot nif Go,"dx f° G* ¢:dx 
—43s —js 
o hs hs 
te 2. ¥ncoth ym J Gom'dx J G* bndx, 
m=2 —tks —}s 


a is 
- — tan Bed r Hdx H*dx 


—js v%—js 


© ts ts 
(27) = my wy HH expli2nx/s)xlde J H*exp[—j(2nm/s)x]dx 
hs 
+6, ‘tan p,l oe ‘ix H* ¢,dx 


+> (ymcoth Yml)~ in Hontds f H* bndx, 
—4s —}s 


where e~/67G (x) is the electric field in the aperture plane z = 0 when the slotted 
section is terminated by a short circuit at z = / and e~/8*H(x)/jwu is the x 
component of the magnetic field in the plane z = 0 under the same conditions. 
Eqs. (26) and (27) are positive definite quadratic forms and therefore yield 
upper and lower bounds on cot Bod respectively. In practice the functions G 
and H are approximated by a finite series. For instance, let 

N 
(28) G = D> -*- exp[—j(2nx/s)x]; 

—N V/s 
then 


N N 2 ~ N 2 
(29) —co Bocot Bod =12,' re +( dX caPu) Bicot Bil+ » ma 3 csPan) , 


where / has been assumed large enough so that coth y,/ may be replaced by 
unity for all m > 1, and where 
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4s }s 
P= = Le omexp|j(2n2/s)x]dx = = £. dm exp|—j(2n2/s)x]dx. 


The general formula for P,,m is given in Appendix II. Since the slotted section 
is terminated in a short circuit all the coefficients c, may be taken as real. 
When the slotted section is terminated in a matched load the coefficients c, 
would, in general, be complex. Furthermore the variational expressions for 
the input admittance and impedance would involve the aperture field for a 
wave incident at an angle —@; in place of the conjugate fields. The coefficients 
cn in Eq. (29) are determined by solving the set of equations obtained by 
setting all the partial derivatives with respect to c, for | | > 0 equal to zero. 
There is no loss in generality by taking cy) equal to unity. The summation 
over m is not terminated but in practice the cross coupling factors P,» approach 
zero rapidly for large m when m # n. For m = n, Pam approaches unity for 
large m. 
Equation (27) may be handled in a similar manner. 


6. NUMERICAL EXAMPLE 
The equivalent circuit parameters for a slotted interface having the fol- 
lowing dimensions: 


S$ = ] cm., 
t = O3d:cm:, 
«= 25e 


were evaluated for a free space wavelength of 3.14 cm. and an angle of incidence 
of 30°. For this example the value of 6; is 2.32 and 8) = 1.73. Sufficient 
accuracy was obtained by using a constant term only for G and HZ. A plot of 
Bil+-Bod vs. Bil is given in Fig. 3 from which it is seen that the maximum error 


2 


B, 4+ Bod 


-d 





“2 = 


1 % 


oO 
Bz-7 
/ 
Fic. 3. Theoretical curve of Bod +,/ vs. 8y/. 


in Bod does not exceed +0.015 radians and hence a position d of a field minimum 
is correct to within +0.009 cm. The turns ratio m of the ideal transformer is 
given by 
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/— 


n= A/ cot ZL™ = 1.004.005, 
Bo 4 
where w is the total peak to peak width of the curve. From the curve the 
electrical lengths of the two transmission lines, « and v, are estimated to be 
0.85°+.4° and —1.15°+.4° respectively. It is seen that there is a negligible 
phase shift introduced into the reflected and transmitted waves. Also the 
impedance transformation at the junction may be taken as that due to the 
change in wave impedance in going from free space into the slotted section. 
The reason for this behavior can be readily appreciated by referring to Fig. 4 


x 





‘oo 4 6 8 


Fic. 4. Comparison of functions F\(x) and e~/8?, 
1—Real part of Fi(x), 2—Imaginary part of F,(x), 
3—Real part of e748, 4—Imaginary part of e~/8?, 


where the real and imaginary parts of F; and e~? are plotted. It is seen that 
F, differs from e~/8? by less than +12% and hence the amplitudes of the 
higher order modes excited are small. 

Since such a slotted dielectric behaves like a homogeneous dielectric with an 
effective dielectric constant x, equal to (6;7/y*)+sin? 6; at the wavelength Apo 
and angle of incidence @;, the design of a slotted dielectric matching layer is 
relatively straightforward. The requirements are that the spacing s must be 


less than 
No 


V/«(1+|sin 6;|) 
and the thickness be such that 
ke—sin’@, = [(1—sin26,) (x—sin’6,)]}. 
The depth of the slots D is chosen as 
Xo 





tnfe,—ain, 
or one quarter wavelength. These equations are the same as for a homogeneous 
dielectric of dielectric constant x,. For the following parameters: 
s=l1cm., t = 0.35 cm., 
D = 0:635 cm., No = S14 emi, 
x = 2.56, 
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a wave incident at an angle of 15° is perfectly matched to a dielectric of relative 
dielectric constant 2.56. The magnitude of the reflection coefficient for angles 


of incidence of 0°, 15°, and 30° in the wavelength range of 3.00 to 3.28 cm. is 


‘05 


‘025 





3:00 3:07 3-14 321 328 
A> IN CM. 


Fic. 5. Modulus of reflection coefficient for a slotted dielectric matching section. 


plotted in Fig. 5. The reflection coefficient does not exceed 0.05 as compared 
to 0.23, 0.24, 0.275 for the above angles of incidence when no matching section 


is used. 
7. CONCLUSIONS 


A variational solution for the parameters of the equivalent circuit of a 
slotted dielectric interface for a wave incident at an angle to the interface 
normal and polarized with the electric vector parallel to the slots has been 
given. For the usual values of dielectric constants and angles of incidence that 
would be encountered in practice the phase shifts introduced into the reflected 
and transmitted waves are very small being of the order of 1°. A quarter wave 
matching section may be constructed from such a slotted dielectric and will 
match a plane wave to the solid dielectric medium over a wide frequency band 
and a considerable range of angles of incidence with negligible reflection. The 
slotted dielectric has some similarity with the parallel plate media. One out- 
standing difference is the possibility of choosing a spacing s such that there 
will only be one reflected wave. For parallel plate media the spacing must be 
greater than a half wavelength so that for angles of incidence above a certain 
critical angle two reflected waves appear (Berz 1951). 

For a wave incident with the magnetic vector parallel to the slots the wave 
impedance and also the wavelength in the slotted section is different so that a 
match will not be obtained. The same is true for a homogeneous dielectric when 
6, ¥ 0 with the exception that the wavelength does not change with a change 
in the polarization. In order to match both polarizations simultaneously an 
additional degree of freedom in the design must be introduced. This can be 
done by cutting an additional set of slots perpendicular to the first set. The 
theoretical analysis of such a structure would be difficult. Some experimental 
work on such a matching section has however been done (Jones and Cohn 
1955). 
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APPENDIX I 
With reference to Fig. 6 the electric field for the m'th mode will be a linear 


combination of two waves propagating in the positive and negative x directions. 
Let the amplitude of the forward and backward travelling waves be a; and }, 


Leese 9 4 


te > : 
We GB 
vA 44 EC 
| 


—_ a 
t RT Ro I 
Fic. 6. Periodic slotted dielectric medium. 


respectively at x = 0 and a» and b» respectively at x = s. The electric field of 
the forward travelling wave is proportional to e~/'™? in the dielectric sections 
and e~/*»* in the free space sections. The z component of the magnetic field is 
proportional to /,e~?'"* and k,e~*"* respectively. For the negative travelling 
wave the signs in the exponential terms are changed to positive. The wave 
admittances may be taken proportional to /,, in the dielectric sections and k» 
in the free space sections. The reflection coefficient R; and transmission co- 
efficient 7; looking to the right at x = 0 are given by the following expressions: 


_ kal 
R m_ = =m 
(30) aT ie 
T; = 14+ R. 
At x = ¢ the corresponding quantities are 
R; = —R, 
(31) 
T> = 1—R,. 


It may be shown (Collin 1955) that the amplitudes a,, b; are related to a2, b2 


by the following matrix equation: 
ay e?: im t+Kma) l Rye J2lmt 1 —R,e a2 
(32) = TT 
b, ile R; en zimt aii, en f2kma b» 
a (4 -— 
An AxJLb: 
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In order to have a mode propagating down the periodic structure in the positive 
x direction it is necessary that a, = yi and b, = e~**b,. Therefore 


e #%q 
by 21 Hy, 
which has a solution if the determinant of the resulting two simultaneous 
equations vanishes. The vanishing of the determinant gives 
(34) A114 22— A 12A 21+ 6748 — €8*(A1,4+A22) = 0. 
When the appropriate values for the matrix elements A ;; are substituted from 


Eq. (32) into Eq. (84) and the resulting equation simplified the following 
eigenvalue equation is obtained: 


2 
(35) cos Bs = cos /,t cos kma — In ool sin J, sin Rm. 
This equation together with the relation ky? = 1,?—(k—1)ko? permits a 
solution for /,, and k, to be obtained when 8 is given. When m becomes large 
lL, and k,, approach equality and the term (/,?+-Rm?)/2/mkm approaches unity. 
The eigenvalue equation thus becomes 

cos Bs = COS(Imt+kma) = cos l_s 

and hence /, = (2mz/s)+8. Under these conditions @,,“° approach zero, a,,°° 
approach unity, and B, approaches jA, for m even and —jA, for m odd. 
Thus the functions ¢,, for m even approach exp[j(2mz/s)x] and for m odd 
exp[—j(2mx/s)x]. When 8 = 0, Eq. (35) reduces to 
(36) [Jn tan I_(3t)+Rm tan Rm(3a)] [lm COt lm (4t)+Rm COt Rm(4a)}] = 


The vanishing of the first bracket gives the eigenvalues for normal mode 
functions which have even symmetry about the points x = 3t, 3(s+42), while 
the vanishing of the second bracket gives the eigenvalues for the normal mode 
functions which have odd symmetry about the points x = }t, 3(s+t). The 


cos(£s sin@;) 
o - 





° 2 4 6 8 10 12 


Fic. 7. Curve for determining the signee l, as a function of angle of incidence, 
t = 0.35 cm., s = 1 cm., ko = 2, « = 2. 
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latter modes are not excited in the case of normal incidence when 8 = 
When 8 # 0 the modes are not separated into two discrete groups. 

Fig. 7 gives a plot of the right-hand side of Eq. (35) as a function of /, for 
s = 1em., t = 0.35 cm., «x = 2.56, and ko = 2. The eigenvalues /,, are readily 
determined from the points of interception of this curve with the straight line 
cos(Ros sin 6;) when the angle of incidence 6; is given. From the figure it is seen 
that the right-hand side of Eq. (35) approaches cos /,s rapidly. 

The eigenvalue equation for a wave polarized with the electric vector in the 
xz-plane may be derived in a similar manner. It is as follows: 


(37) cos Bs = cos 1,,t cos ka * fe le. sin Jjt sin Rma. 





This equation does not have a simultaneous solution with Eq. (35) so the 
wavelength in the slotted section will be different for the two polarizations. 
When 6 = 0 the normal modes again split into two sets having even and odd 
symmetry properties and with the following eigenvalue equations respectively : 


(38) ln tan ly (3t)+xkm tan km(4a) = 0 
Im COt Im (3t)+kRm COt Rm (Fa) = 
APPENDIX II 
The coefficients A,, are chosen so that the normal mode functions F,, are 
normalized. Thus 
4s 
(39) FF gtd = An 7. L(Y’)? ten? (Ym?)* Hex = 
—}s 
Substituting for ¥,,° and y,,” and integrating gives the following expression 


for Am 


~ 


An = = (1+am° 43 [(am’)” +0m" (am’)’| 


2 





(40) 1 Gm > isin Qkm(45—Om°) —sin 2km(4t—Om°)] 





— an") [sin 2km(#5—On”) —sin 2m (4 —Om”) |. 


The cross coupling factors Pym are defined by the following integral 


hs 
Pa J dm exp[j (2am /s)x|dx 
-is VS 


‘ hs ‘ 
— J | va‘cos( 04-28) tanibn’sin( +222) lee 


(41) 
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Carrying out the integrations gives the following expression for Pym: 
1+a, 
ie 


n 





(42) PamvS D8 ain rng (Ht) +E sin Hina (3) 


mn 





+a + sink (3s) oa mm | — sin [foam ( 3t) ae RmPrn’ | 








a 

4 Sin[gnn (38) — Ridin’ — sin [gna (41) bata! | 

&mn 
= kn Om’) — sin[gmn(3t) — kmOn’] 
Bmn 

__sit[fran (45) — kmOm?]—sinn [fran (32) — bata’) 

Sun . 
where mn = lm+B+2nz/s, 


inn = In —B8—2nx/s, 
Sun = km+B+2nx/s, 
Lnin = Rm—B—2nw/s. 





THE AXIAL-FLOW ERROR IN THE THERMAL- 
CONDUCTIVITY PROBE! 


By J. H. BLACKWELL’ 


ABSTRACT 


In connection with the development of the thermal-conductivity probe the 
author has used a new boundary condition of Jaeger to find an expression for an 
upper limit to the axial-flow error in these devices. Unlike in previous attempts 
at this problem the short-circuiting effect of the probe material is taken into 
account. 


INTRODUCTION 


The idealization of “radial flow” is frequently employed in the solution of 
transient heat flow problems with axial symmetry in an infinite medium but 
convenient criteria have not often been available for determining the minimum 
length of the source or sink of heat for this idealization to be valid. 

A recent paper by Jaeger (1955) makes possible solutions of this problem 
for the important practical cases where the heat source is a thin wire or thin- 
walled hollow cylinder of good conductor. 

It is the purpose of the present note to show how Jaeger’s radial boundary 
condition may be used to find, in convenient form, an upper limit to the axial- 
flow error in thermal-conductivity probes (Blackwell 1954; Hooper and Lepper 
1950). 


DESCRIPTION OF THE PROBLEM 


A thermal-conductivity probe is normally a circular cylinder of good thermal 
conductor. It is either solid with a small radius or hollow with a thin wall so 
that radial temperature differences in the probe material are negligible to 
first order. It is provided with a heater and some means of measuring tempera- 
ture at the center of its length. 

In practical applications the probe is either buried in the external medium 
whose thermal constants are being measured (granular material) or inserted 
in a long hole (solid material). In the design of these devices it is obviously 
important to determine in advance the minimum length of the probe for 
simple radial flow theory to be valid. 

It is possible to set up idealized radial—axial heat flow problems which must 
have greater or smaller axial flow effects than occur with a practical probe. 
Case II of a previous paper by the author (1953) provides an example of the 
latter and it is the present purpose to solve a corresponding example of the 
former. 

The problem considered is that of an infinitely long cylinder of good con- 
ductor, whose cross-section satisfies the conditions given above for a practical 
~ 1Manusc ript 1 rec eived November 30, 1955. 

Contribution from the Department of Geophysics, Research School of Physical Sciences, 


Australian National University, Canberra, Australia. 
*Visiting Research Fellow; on leave from University of Western Ontario. 
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probe, and which is immersed in an external medium of infinite extent. Heat 
is supplied to a finite length of this cylinder only, corresponding to the length 
of the probe. 

Then, provided that a real probe is immersed sufficiently well that external 
boundary effects are negligible and that other mechanisms of heat transfer 
than conduction may be disregarded, the real axial heat loss must be less than 
that of the idealized problem. With the possible exception of the presence of 
convection in certain geophysical applications, the restrictions stated have 
little experimental significance. 


ANALYSIS 
We wish to solve the following problem and compare solutions with the 
corresponding one for radial flow: 


(1) oot 90, 98 _ 1 a<p<m~, —x<zs<om, t>0, 
p Op dz’ «OF 


subject to 


(2) ¢é=0, t=0,a<p<o,—-~7 <2z< oO, 
6, = 0, t= 0; 


(3) -K5 = aK A to7¢64-2)—- ue 1) + 54-4 21] 


= H(6,—6), p=a,t>0, 


where 
p,2,t are the radial and axial co-ordinates and the time, 
6(p, z,t), K,x are the temperature, conductivity, and diffusivity of the 
external medium, 
6,(z, t), Ki, x, are the corresponding quantities for the probe, 
a__is the probe radius, 
L _ is the probe half-length, 
5 _is the probe cross-section, 
A _ is the heat supplied to the probe/unit volume/unit time, 
H is the “outer conductivity” at p = a (a measure of the 
contact resistance between probe and medium), 
U(x) is the unit step-function. 
Equation (3) is Jaeger’s boundary condition (equation (21) of his paper) 
modified for the presence of contact resistance and to include the axial con- 


ditions of the present problem. 
By standard methods, we find the Laplace transform of the temperature 


in the cylinder: 


bay aah 1) ly sin Lw e'™{Ko(8a)+(K /aH) Ba K;(8a)} dw 
pRy 2a 


~« «| (6'a°+6) {Ko(8a)+ (K/aH) Ba K,(8a)}+BfaK,(8a)) 
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where pis the transform variable, 
= (q? + w)}, 
q = (p/«)}, 
b = wa?(m/x—1), 
B = 2na? Kni/6Kix. 

The axial-flow error obviously increases with time and in the case of pure 
radial-flow (Blackwell 1954) the conductivity K is inversely proportional to 
06;/A(In t) for T, = xt/a* > 1. It is therefore permissible and also convenient 
to discuss the axial flow error in terms of an approximate value of 00,/0(In ¢)],=0 
valid for large 7,. If in our analysis, we take into account the further require- 
ment that L/a is large, we arrive at a simple relationship between 06,/d(In ¢) 


in the radial and radial—axial cases. 
By standard methods of the Laplace Transformation and assuming orders 


of integration may be reversed, we obtain from equation (4), 


© LL = Be) SL Sento Ferd 


where 
(6) F(w) = {Ko(qa)+(K /aH)qaK (qa) } ed 
ae pr: (qa +b) {Ko(qa)+(K /aH)gaK 1 (qa) } +BqaK (qa) 


and Br2 is the contour ABCDE of Fig. 1. 


Fic. 1. 


Following Goldstein's method we then expand the integrand of F(w) in 
ascending powers of ga; term-by-term integration should then result in an 
expansion suitable for large 7, (Blackwell 1954; Goldstein 1932). In the present 
instance we retain dominant terms only and get 


F(w) = _{In(}qa)+7—(K/aH)} e"dp_ 


(7) = arr —B+b{In(4qa)+y—(K/aH)} 


where y is Euler’s constant = 0.5772. 
Since the major contribution to the integral of equation (5) comes from 
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values of w such that | wl | « =z, we can treat wa? and hence b/B as small 
quantities, and expand the integrand of (7) again, obtaining 


(8) Fw) = ( -5 Bea [ma 29a) +y— K4(2)iing 3qga)+y— x 


+( a) atta | e"dp 

af aH) *°"° 
The term integrals of equation (8) are either well-known or easily derivable 

(Blackwell 1954; Jeffreys and Jeffreys 1946) and we find 


1 b 2K 
0) Bh ($)fnatnr] 


+f 2 —tvalnan— ee teen 


Introduce the dimensionless parameters 
€ = K,/K, 


n = Kyx/Kn, (ratio of thermal capacities/unit volume), 


and o = 42) = 2d for hollow probes, 
a 2a a 











=1 for ‘‘needle’’ probes. 






Then 






6 = ma’o, 
B = 2/on, 
b/B = 3u*a*(e—n), 
and substituting en (9) into Or 


2 © 
‘ pha <-9) w sin Lw exp(—w' xt) dw 


2a Bs 


2 42 2 ro) 
+{ ~ 2 4(in 4T ey jeteest f w'sin Lw exp(—w' xt) ao] ‘ 


00; = 2.) 
(11) ioe |, = (& R 


= erfx —2-2CSxte-2* + 2 -4C2{ — hn? 4357} x4(4+ 4x —x2)e-*?; 
L/2V/xt, C = oa*(e—n)/L’, 
= In 47,-—y+2K/aH. 






















& 
" 







The factor R gives the ratio between the respective values of 06,/d(In #) 
and hence the required limit to the axial-flow error. The first term of R is 
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independent of the probe parameters and may be thought of as due to axial 
flow in the external medium alone; for a given value of ¢ it approaches unity 
as L — o. The remaining terms in the bracket correct the first for the presence 
of probe and contact resistance and approach zero as L — . For design 
purposes, we need only retain the second term, the third giving an estimate 
of the error made in assuming b/B <1. As expected the second term varies 
directly with the ratio of conductivities and the probe cross-section. 

As 6 = ma*e > 0, R—-erfx, a result which happens to be true exactly 
(Blackwell 1953). This value of R provides the lower limit to the axial-flow 
error referred to earlier. 

The parameters C and x used above both contain L. It is convenient for 
design purposes to introduce a further parameter 


\ = L/a, the length-diameter ratio. 


Rewriting the first two terms of equation (12) 
(13) R = erf[A(4T,)-4] —20n-*A(e—n) (47,3? S exp[ —A2/4T,]. 


Inserting the asymptotic expression for the error function we obtain the 
relative error due to axial flow 


(14) AR = 1—R = w+ exp[—*/4T,] [(47,)4/A+20d(e—n) (4T+)-*75] . 


The error in AR introduced by this further approximation is <0.0004 for 

r(4T,)-? 22. 
NUMERICAL EXAMPLE 

The maximum value of 7, used in an experiment is set by consideration of 
the basic radial-flow theory and is unlikely ever to exceed 25. If we take this 
value for T,, and 2K/aH = 2 (a rather large contact-resistance for many 
practical applications), equation (14) becomes 
(15) (AR) max = exp[—0.01A?] [5.64/A+0.680X 10-%cA (e —7)] . 


Equation (15) is an appropriate form in which to use the criterion. 

Take the case of a hollow brass probe of, say, 1} in. O.D. X ¢ in. wall for 
use in a geophysical application where we know K ~ 0.005 c.g.s., «~~ 0.01 c.g.s. 

How long should the probe be? 

We have e = 50.0, 7 = 1.50, o = 0.40, and hence 


(16) (AR) max = [5.64/A+0.132d] exp[ —0.01A2] . 


The expression (16) is now evaluated for trial values of X: 


» 5.64/d 0.132d exp[—0.01A4| (QR) max 


20 0.282 : 0.0183 0.053 
25 0.226 3.30 0.00192 0.007 
30 0.188 3.96 0.000123 0.00051 


We observe that \ = 25 gives a maximum error slightly under 1% while 
for \ = 30, the error is negligible for any normal thermal experiment. 
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On the other hand, consider a solid probe of the same outside diameter. 
The parameter o is now unity and for 


d = 25, (AR) mas = 0.017; 
A= 30; CAR) aa = GOOl2. 
It will be noted that the relative errors have been more than doubled. 


ACKNOWLEDGMENTS 


The author wishes to express his gratitude to Prof. J. C. Jaeger for helpful 
comments throughout. He would also like to thank the referee for suggesting 
the final approximation contained in equation (14) above. 


REFERENCES 


BLACKWELL, J. H. 1953. Can. J. Phys. 31: 472. 
1954. J. Appl. Phys. 25(No. 2): 137. 
Go.psTEIN, S. 1932. Proc. London Math. Soc. 34: 51. 
Hooper, F. C. and Lepper, F. R. 1950. J. Am. Soc. Heating Ventilating Engrs. 22: 129. 
i J.C. 1955. Quart. J. Mech. Appl. Math. 8(Pt. 1): 101. 
EFFREYS, H. and JEFFREYS, B. 1946. Mathematical physics. Cambridge University Press, 
London. p. 373. 











NOTES 


INTEGRALS OF INTEREST IN METALLIC CONDUCTIVITY* 
By D. K. C. MAcDOoNALD AND Lois T. TOWLE 


The electrical resistivity, p, of an idealized metal may be expressed in the 


form: 
p re ink oad 2*dz 
(1) St ee rf nie) 
Po Ta OSe (€ —1)(1-e%) 


. 
: e 6° 


an 


_ J5(0/T), 


a formula due originally to Bloch (1929, 1930). p,, is the resistivity at some 
sufficiently high temperature 7., (>6), and 6 is the characteristic (Debye) 
temperature of the metal. 

The corresponding expression for the electronic thermal resistivity, W, is 
considerably more complex (cf. Wilson 1953, p. 285) and is in poor agreement 
with the experimental data (cf. e.g. Wilson Joc. cit., pp. 288-290). A crude 
physical argument suggests that corresponding to equation (1) we might 
expect the thermal resistivity to be given approximately by: 


WwW 7° ow z*dz 
9 _—_—_- = —_—_————— 
(2) W., 2B J, (e?—1)(1—e*) 


2 
= 24, J:0/T) 


and indeed White, Woods, and MacDonald (1956) have found that this ex- 
pression fits the experimental results surprisingly well. 

Sondheimer (1950) has tabulated the integrals J,(x) with r = 5, 7, 9, 11, 
13, 15, 17 for x ranging from 0.8 to , and Griineisen (1933) earlier gave a very 
full table of (4/x*)J;(x), this form being chosen to have the limit unity as x 
tends to zero. It appears possible that values of J2(x) and J.(x) may be useful 
in the future as well as J3(x) which occurs in (2) above, and we have therefore 
tabulated here these integrals for x from 0.1 to » (J,(x) does not converge at 
the lower limit). For completeness we have also included J4(x) which arises in 
the Debye formula for the specific heat C, of a solid: 


rrr z‘dz 
© c= 9R( 5) Jl, waa 


won 2) 2( 4) 


*Issued as N.R.C. No. 3895. 
{We do not here enter into any discussion of the appropriate choice for 6. 
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Wain detdhtenes (0 cas 
TABLE OF J,(x) ~f (—na—e~) 








x Jo(x) J3(x) J a(x) J5(x) 
6.1 0.09997 0.0050 0.0003 0.000002 
0.25 0.2496 0.0312 0.0052 0.00019 
0.5 0.4966 0.1237 0.04115 0.00616 
1.0 0.9730 0.4798 0.3172 0.1885 
1.2 1.1540 0.6788 0.5365 0.4573 
1.5 1.4122 1.0269 1.0079 1.3319 
2 1.8017 1.706 2.2016 5.0858 
3 2.4111 3.211 5.9632 29.69 
4 2.8067 4.579 10.7293 88.59 
5 3.0392 5.614 15.3671 182.22 
6 3.1657 6.3033 19.1210 295.40 
8 3.2624 6.9581 23.5874 507.00 

10 3.2844 7.1505 25.2812 639.21 
13 3.2895 7.2061 25.9273 713.85 
20 3.2899 7.2123 25.9639 732.31 
@ 


3.2899 7.2124 25.9757 732.49 


: NoTE: Lim J,(x) = x™'/(r—1); - 
z—0 
Lim J,(x) =r! z 4 =r! ¢(r), where ¢(r) is the Riemann zeta-function. 
© s=1 


Hence J2(o) = x?/3, 
Ji(~) = 44/15, 
Je6(~o) = 16r*/21. 
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ON THE THERMOELECTRIC POWER OF DILUTE METALLIC ALLOYS* 
By A. B. Buatrat AND D. K. C. MACDONALD 


The thermoelectric power S of a metal is given by (e.g. Mott and Jones 1936) 
oe vET) Og aE) 


(1) ieee aE Seu 
o ~#ETS ab o(E)\ 
iv 3et dlog E Jee’ 


where p(£) is the electrical resistivity corresponding to some value E of the 
Fermi energy, ¢ is the actual Fermi energy, and the other symbols in (1) have 





*Issued as N.R.C. No. 3909. 
fNational Research Laboratories Postdoctorate Fellow; now at Physics Department, 


University of Alberta, Edmonton, Alberta. 











420 CANADIAN JOURNAL OF PHYSICS. VOL. 34, 1956 


their usual significance. If then we assume Matthiessen’s rule to be valid we 
may write: 
(2) P = prt po, 


Caen Sean 


where pr is the resistivity of a metal in the absence of impurities, due to the 
thermal agitation of the metal ions, and pp is the resistivity due to chemical 
impurities and physical imperfections. Using (2) in (1), we then have: 
(3) oo eee _ 9 log p Pr, Be _ 9 log a | 
3ee Lp Olog Ep @log Edguy 

Further discussion in this note will be confined to the quantity dlogp)/dlogE 
which for brevity we call C. 

If we assume quasi-free electrons (E « | k |?), then 


ara ine A) 
(4) or 1+(2 log E/ g.;’ 


where A is the effective scattering cross-section per impurity center. In the 
case of an impurity atom of valence Z + 1 dissolved in a monovalent metal, 
Mott (1936) calculated the potential approximately* due to this atom using 
the Thomas-Fermi method, and obtained: 


(5) V(r) = (Ze/r) exp(—r/ro), 
where 7 is a screening radius which in this treatment is independent of Z. 
Using the Born approximation, Mott then found: 

; — Pe E/E, \ 
(6a) A= oR log(1+E/E1)—TT 7B, : 
where E = mey,v?/2, Ey = h?/8mertro?, Mere being the effective electron mass. 
If now we treat E; as a constant for the metal we obtain for C the curve of 
Fig. 1 as a function of ¢/F;. It will be seen that C lies between —1 and —3 





Fic. 1. Variation of C using screened Coulomb potential. 


*Fujiwara (1955) has recently obtained a numerical solution of the exact Thomas-Fermi 
equation. 
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(cf. MacDonald and Pearson 1953, and also Domenicali and Otter 1954). 
Alternatively we may substitute for ro the theoretical value dependent on 
electron-concentration derived by Mott, and we have (still assuming quasi- 
free electrons): 


_ 7Zie'f > /m\9..._(E/Es)* \ 





where Eo = mMezre!/27°h?, ~1.35 ev. if mere ~ m. The variation of C in this 
case is also shown in Fig. 1 as a function of [/E2, and C now lies between —2 
and —3. 

However, measurements in this laboratory (e.g. MacDonald and Pearson 
1953, 1954, 1955) of the thermoelectric power at low temperatures of certain 
low concentration alloys (e.g. tin and iron in copper) indicate a much wider 
range of values of C (sometimes |C| may even exceed 100) than obtained from 
Mott’s formula. If we assume the validity* of equations (1) to (3) and of the 
assumption of quasi-free electrons, the explanation for this discrepancy must 
be sought in the calculations of V(r) and of the scattering cross-section A. 

Since the screened Coulomb potential (5) for an impurity ion is only a rather 
crude approximation which in particular takes no account of the specific ion 
structure, we thought it worth while to consider also the potential V = 1, + V2 
where: 


(7) Vi = (A/r) exp(—r/ro); V2 = BforO cr gn, B= —(8A/mn)(r0/n1)?, 
= O'forr > r:. 


The calculation of scattering cross-section is straightforward on the Born 
approximation, and the range of variation of C is greatest if we set ro/r1 = 1/38. 
C is plotted against kr; in Fig. 2 (where & is the electron wave-number) and 
we see that C now lies between +3 and —3. If we assume r to be about the 
radius of the Wigner-—Seitz atomic sphere, the value of kr; at the Fermi surface 
will be about 1.5 and hence C ~ +1. Thus although the possible range of C 
has been increased, we are still as far as ever from obtaining values of C com- 
parable with those observed experimentally. 

We should also mention that Friedel (1955)f has shown that when the 
potential V of the impurity center is such that there exists a virtual state(s) of 
some non-zero angular momentum close to the Fermi energy, such a state(s) 
can cause the appearance of a maximum and minimum in the p(£) vs. — curve. 
Somewhat higher values of |C| would now be possible and Friedel has shown, 
assuming a square-well form of potential, that the existence of such a virtual 
state is not incompatible with a reasonable choice of parameters for the well. 
It still seems improbable, however, that we can increase C by some two orders 


*It appears that any considerable modification of the E-k relation from that assumed can 
be ruled out as an adequate explanation: (i) since this would effect a large change in the 
thermoelectric power of the pure parent metal also (cf. e.g. MacDonald and Roy, Phil. Mag. 
44: 1364, 1953), and (ii) because experimentally the values of |C! depend markedly on the 
specific solute atom. 

‘ tWe are most grateful to Dr. Friedel for sending us a copy of his manuscript before pub- 
ication. 
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+3 
+2 


+1 


kr, 
Fic. 2. Variation of C using compound potential specified in text. 


of magnitude. An exact calculation to determine the upper bound of |C! for a 
square well potential is in progress. 

We should emphasize in conclusion that any explanation based on equation 
(2) assuming a linear variation of pp with impurity concentration can only be 
incomplete, particularly at low temperatures. This is particularly so because 
of the striking ‘‘saturation” effects observed (e.g. MacDonald and Pearson, /oc. 
cit.) in certain very dilute alloys which appear undoubtedly to be intimately 
connected with the occurrence of an anomalous minimum of electrical re- 
sistance in these alloys at low temperatures (cf. also MacDonald 1953, 1955). 
We therefore believe that a theory should continue to be sought for a specifi- 
cally low-temperature component of electron scattering which could give rise 
to the anomalous resistance-minimum and large thermoelectric power in 
dilute alloys.* 

We are grateful to Dr. W. B. Pearson, D.F.C., for discussions. 
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*We should like to mention here the work of Korringa and Gerritsen (1953). Let us remark 
also that recent work by Pearson (1955) has enabled the anomalous component of resistance 
to be isolated experimentally and this shows indeed a rapid decay with increasing temperature 
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